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Abstract 

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity 
of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar 
mass range of the parent YSOs. To this end we have obtained optical to near-infrared spectra of HH 1042 and HH 1043, two newly 
discovered jets in the massive star-forming region RCW 36, using X-shooter on the ESO Very Large Telescope. HH 1042 is associated 
with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra of both targets. High-velocity (up to 
220 km s~^) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low- velocity emission from the 
background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared 
spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass 
outflow rates are Mjet ~ 10~^ M© yr~^ It is not possible to determine a reliable estimate for the accretion rate of the driving source 
of HH 1043 using optical tracers. We measure a high accretion rate for the driving source of HH 1042 (Mace ~ 10~^ M© yr"^. For 
this system the ratio Mjet /Mace ~ 0.1, which is comparable to low-mass sources and consistent with models for magneto-centrifugal 
jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. 
While the mean velocities in both lobes of the jets are comparable, the variations in mass outflow rate and velocity in the two lobes 
are not symmetric. This asymmetry suggests that the launching mechanism on either side of the accretion disk is not synchronized. 
For the HH 1042 jet, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description 
of a ballistic flow as its constituents. We have simulated the flow and the resulting emission in position- velocity space, which is then 
compared to the observed kinematic structure. The knotted structure and velocity spread can be reproduced qualitatively with the 
model. The results of the simulation indicate that the outflow velocity varies on timescales on the order of 100 yr. 

Key words. Stars: formation - Stars: circumstellar matter, interstellar medium (ISM) - ISM: jets and outflows - ISM: Herbig-Haro 
objects - ISM: individual objects: HH 1042 - ISM: individual objects: HH 1043 



1. Introduction 

Astrophysical jets are a ubiquitous signature of accretion. When 
the magnetic field of the circumstellar disk is coupled to a 
jet, magnetic torques can remove a significant fraction of the 
system's angular momentum and mass through the jet, when 
charged particles are ejected. Jets exist in accreting systems on 
various scales, ranging from young stellar objects (YSOs) up to 
evolved binary systems and active galactic nuclei (AGNs). They 
are detected in emission lines at all wavelengths, from the ra- 
dio domain up to X-rays. Jets associated with optical emission 
are classified as Herbig-Haro (HH) objects. The velocities mea- 
sured in jets are close to the escape velocity at the launch region: 
from a few hundred km s~^ in YSO jets up to relativistic veloci- 
ties in AGN jets. The mass flux also scales with the mass of the 
central object (Miley, 1980; Bally et al., 2007). 



* Based on observations performed with X-shooter (program P87.C- 
0442) mounted on the ESO Very Large Telescope on Cerro Paranal, 
Chile 



Most jets show varying velocities and shock fronts along the 
flow axis, which can be attributed to a launching mechanism at 
the jet base that is variable in time (Rees, 1978; Raga et al., 
1990). As suggested by Reipurth & Aspin (1997), the shocked 
structure of HH jets may reflect an FU Orionis-like accretion 
process: relatively short periods of intense accretion, which may 
be caused by thermal instabilities in the accretion disk. In prin- 
ciple, then, a history of the accretion is contained in the fossil 
record of the shock fronts in the jet. Jets thus provide the oppor- 
tunity to obtain time-resolved information on the variability of 
the accretion process from a single observation. Given the typ- 
ical spatial extent (up to a few pc) and velocities of jets it is 
possible to probe accretion variability on dynamical timescales 
up to a few thousand years, much longer than what is possible 
with a series of direct observations of the accretion process. 

The jet launching mechanism described in the seminal paper 
by Blandford & Payne (1982) serves as the standard model of 
jet formation. Matter is removed from the accretion disk by cen- 
trifugal forces, and confined along magnetic field lines that carry 
away the material from the disk in bipolar directions. Pudritz & 
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Norman (1983) proposed that the same mechanism drives YSO 
jets, as the disks of young stars are known to be threaded by mag- 
netic fields. Most of the current understanding of proto stellar jets 
comes from studies of low-mass objects, since these are much 
more numerous, and they form on longer timescales than their 
higher mass counterparts. However, collimated outflows are also 
detected around some massive YSOs (e.g. Torrefies et al., 201 1). 
As the luminosity of the central star increases, radiation pressure 
is thought to become a more significant driver of jets (Vaidya 
et al., 2011), in addition to the centrifugally driven magneto- 
hydrodynamic disk wind. 

Just over a dozen forming (potentially) massive stars (M > 
10 Mo) have been reported, most of them through the detec- 
tion of disk signatures at infrared, sub-mm, or radio wavelengths 
(e.g. Cesaroni et al., 2006; Zapata et al., 2009; Davies et al., 
2010; Torrefies et al., 2011). The observed rotating disks around 
these objects can often be associated with collimated molecu- 
lar outfiows. Jets from massive stars have been observed for a 
few objects in mm tracers, with estimated mass outfiow rates in 
the range of 10"^ to 10"^ M© yr"^ (Cesaroni et al, 2007). In 
contrast, jets from low-mass stars have been observed mostly at 
optical to near-infrared wavelengths, with mass outfiow rates in 
the range of 10"^^ to 10"^ M© yr'^ (Hartigan et al, 1995; Cofi'ey 
et al., 2008). The highest rates are measured in the less evolved 
sources, i.e. Class I objects (Hartigan et al., 1994; Bacciotti 
& Eislofi'el, 1999; Podio et al., 2006) and in FU Orionis ob- 
jects (Calvet, 1998). The lowest rates are measured in classical 
T Tauri stars (CTTS), i.e. Class II objects. 

There currently have only been a few observations of jets 
from intermediate mass (2-10 M©) YSOs in the optical to near- 
infrared. These are mostly evolved objects, i.e. Herbig Ae/Be 
stars (HAeBe; e.g. Levreault, 1984; Nisini et al., 1995; Wassefi 
et al., 2006; Melnikov et al., 2008), with mass outfiow rates of 
10"^ - 10"^ Mo yr-\ comparable to those of CTTS. As more 
massive YSOs are more embedded, these are the most massive 
objects whose jets can be studied in the optical and near-infrared. 
Their observation may lead to constraints on the jet launching 
mechanism as it scales up to higher masses. 

In this paper we study two jets which were recently dis- 
covered by Ellerbroek et al. (2011). These authors report the 
discovery of a disk-jet system around the intermediate-mass 
(M, ~ 2-5 Mo) YSO 08576nr292, located in the young massive 
star-forming region RCW 36 (Bik et al, 2005, 2006, Ellerbroek 
et al., in prep.). This region is located in the Vela molecular ridge 
at an estimated distance of 0.7 kpc (Liseau et al., 1992). The 
YSO lies at the periphery of the star-forming region, about V 
(~ 0.2 pc) to the West of the center of the cluster, which contains 
two late 0-type stars. The spectrum of 08576nr292 is dominated 
by continuum emission from an accretion disk, with many emis- 
sion lines originating in the disk, in the accretion columns and 
in the outfiow. The jet is spatially resolved by VLT/SINFONI 
integral-field H- and ^-band spectroscopy and shows a clumped 
velocity structure in [Fe ii] and H i emission lines. The veloc- 
ity of these lines coincides with the outfiow indicators in the 
VLT/X-shooter spectrum, suggesting that the jet originates from 
a spatially unresolved region close to the star. The width of the 
jets is not spatially resolved. 

In the same study, the discovery of another jet system is re- 
ported, which emerges from object 08576nr480. Follow-up ob- 
servations of both jets and their sources were carried out with X- 
shooter. The broad spectral coverage and intermediate spectral 
resolution of this instrument provide the opportunity to study the 
physical and kinematic properties of the jet simultaneously in the 
optical and in the near-infrared. The detection of optical emis- 




Figure 1. Detail of the [Fe ii] line map of RCW 36 (d = 0.7 kpc) ob- 
tained with SINFONI (Ellerbroek et al. 2011). Merged slit positions 
during the X-shooter observations are indicated, as well as the positions 
of the knots defined in Fig. 3. Two O stars are located in the central clus- 
ter region, ~ 0.5' eastward. 



Table 1. List of the X-shooter observations. 



Object 


HH 1042/ 


HH 1043 / 




08576nr292 


08576nr480 


a (J2000) 


08:59:21.67 


08:59:23.65 


6 (J2000) 


-43:45:31.05 


-43:45:30.51 


Date 


19-01-2011 


12-02-2011 


Exposure time"^ (s) 


600 


600 


Position angle (°, N to E) 


129 


97 


Sky frame off'set Aa, A6 (") 


44,54 


30,248 


Slit uvb/vis/nir 


1.0/0.9/0.4 


1.0/0.9/0.6 


Resolution Av (km s~^ ) 


59 / 34 / 26 


59/34/37 


^-band seeing (") 


0.6 


0.8 



^ Per slit position 



sion in the jets resulted in their inclusion in the updated version 
of the catalog of Herbig Haro objects (Reipurth, private commu- 
nication) under the entries HH 1042 (the 08576nr292 jet) and 
HH 1043 (the 08576nr480 jet). In this study we use these HH 
numbers when referring to the jets, and maintain the original 
nomenclature (08576nr292, 08576nr480) when referring to the 
central sources. 

The analysis presented in this paper consists of three parts: 
a presentation of the optical to near-infrared spectra of HH 1042 
and HH 1043, a description of the physical conditions in these 
jets, and a simulation of the kinematics of HH 1042. In Sect. 2 
we describe the observations and data reduction, and in Sect. 3 
we present the obtained spectra. Physical conditions and the 
mass outflow rate in the jets are estimated by applying spec- 
tral diagnostics, while emission lines excited in the accretion 
columns provide an estimate of the accretion rate (Sect. 4). The 
kinematic structure of the HH 1042 jet is simulated with an inter- 
pretative physical model; this is presented in Sect. 5. In Sect. 6 
we discuss the appearance of the jets, the accretion and mass 
loss rates, and constraints on the launching mechanism. Finally, 
Sect. 7 contains a summary of this work. 
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2. Observations and data reduction 

X-shooter is a cross-dispersed echelle spectrograph mounted on 
UT2 of the ESO Wery Large Telescope, which produces a spec- 
trum at every spatial pixel along its 11" long slit, in three sep- 
arate arms: UVB (300-590 nm), VIS (550-1020 nm) and MR 
(1000-2480 nm). The slit width can be chosen individually in 
each spectrograph arm (D'Odorico et al., 2006; Vernet et al., 
2011). 

The X-shooter slit was aligned with the jets, in two offset 
positions, both including the star, with a relative offset of 8'' in 
the slit direction, so that the central around the source was 
covered twice (Fig. 1). In between, an exposure was taken on an 
empty part of the sky northeast of the target (Aor = +54'', = 
+44'' NE of 08576nr292; Aa = -30", A^ = -248" SW of 
08576nr480). The obtained two-dimensional spectrum covers a 
wavelength range from 300 to 2500 nm and ~ 19" along the jet, 
i.e. the first 9.5" (~ 6,500 AU at a distance of 0.7 kpc) of the ap- 
proaching ('blue') and receding ('red') lobe. In some cases the 
jet extends beyond the full length of the slit. Table 1 lists the 
characteristics of the observations. 

The raw frames were reduced using the X-shooter pipeline 
(version 1.3.7, Modigliani et al., 2010), employing the stan- 
dard steps of data reduction, i.e. order extraction, flat fielding, 
wavelength calibration and sky subtraction, to produce two- 
dimensional spectra. Observations of the telluric standard star 
HD80055 (AOV) and the spectrophotometric standard star GD71 
(a DA white dwarf) on the same night were used for the re- 
moval of telluric absorption lines and flux-calibration. The ab- 
solute flux calibration agrees to within 3 - 10% with the existing 
photometry. A scaling in the relative flux calibration between the 
two lobes was performed on the NIR spectrum of the HH 1042 
jet, where the use of a narrow slit caused some relative slitlosses 
between the observations of the two lobes. 

The wavelength calibration in the VIS and NIR arms was re- 
fined by using the telluric OH emission lines. The UVB arm was 
subsequently calibrated on the VIS arm by the use of the Na i D 
feature at 589 nm, which appears in both arms. The wavelength 
array was then calibrated with respect to the local standard of 
rest (LSR). 

As shown in EUerbroek et al. (2011) and Fig. 1, the objects 
are embedded in a star-forming region in which the ionizing stars 
have ionized part of the ambient cloud. Throughout the paper, we 
refer to this ambient interstellar medium as the 'cloud'. Since 
photospheric features are lacking in the spectra of the central 
sources due to strong veiling and/or extinction, we assume the 
systemic velocity to be equal to the cloud velocity as measured 
from the nebular lines detected in our spectra (see Sect. 3.1). The 
nebular lines are at -6.5+2.8 km s"^ and -1.0+6.2 km s"^ (^lsr) 
for HH 1042 and HH 1043, respectively. We correct all the spec- 
tra for these values. Thus, the velocities mentioned throughout 
the paper and in the plots are systemic velocities (^sys)? i-^- those 
with respect to the cloud velocity. 

Figs. 2, 3, and 5 show position-velocity diagrams of the 
jets in a number of emission lines. Fig. 4 displays the one- 
dimensional on- source spectra. 

3. Analysis of the emission line spectra 

The spectra obtained for the two HH objects contain more than 
90 emission lines of atomic (neutral and ionized) and molecu- 
lar species at diflferent velocities. They trace various phenomena 
and physical conditions within the system. As we show in the 
following subsections the detected emission lines are originat- 



ing from: (i) the ambient cloud, (ii) the circumstellar disk and 
the accretion/ejection region (on-source) and (iii) the jet. 

3. t. Emission from tine ambient cloud 

After correcting for the systemic velocity (i.e. the cloud veloc- 
ity) some of the lines detected in the spectra have emission cen- 
tered at i^sys = km s"^ (see Fig. 2). This emission appears in 
many lines (e.g., H, C, N, O and S) along the entire slit with 
no velocity variation, suggesting that it originates in the ambient 
cloud. The cloud is not emitting in the ions of refractory species 
that are strongly depleted in the ISM (e.g. Fe n, Ni ii, Sembach 
& Savage 1996). On the other hand, we detect cloud emission 
from highly ionized species (Eion > 40 keV), such as [O m] and 
[Ar III], due to the illumination by recently formed massive stars 
(cf. EUerbroek et al., in prep.) in the central region of RCW 36. 

The emission is also present along the slit in the observations 
of EUerbroek et al. (201 1), where the slit was placed perpendic- 
ular to the jet (see their Fig. 2). The extent of the cloud emission 
is confirmed by the SINFONI Bry and H2 linemaps of RCW 36. 
In Tab. 2 we list the lines that are detected in the cloud (i.e. at 
i;sys = km s~^), along with the lines detected in the jet (i.e., at 
high blue- and red-shifted velocities, see Sect. 3.3). 

3.2. On-source emission: disk, accretion and outflow tracers 

Fig. 4 displays sections of the spectra of 08576nr292 and 
08576nr480 (the central sources of HH 1042 and HH 1043, re- 
spectively), extracted from the two-dimensional frames. Neither 
spectrum contains photospheric absorption lines. The spectral 
energy distribution (SED) of 08576nr292 is dominated by emis- 
sion from a circumstellar disk (EUerbroek et al., 2011). Its spec- 
trum is very rich in emission lines, which trace the circumstel- 
lar disk, a stellar or disk wind, (possibly) magnetospheric accre- 
tion and the onset of the jet. Hot and the Ca 11 triplet lines show 
blueshifted absorption by the jet or disk wind. Although no di- 
rect signature of infall (red- shifted absorption) is detected, var- 
ious emission lines associated with accretion activity (e.g. H i, 
Ca II, and He i) are used in Sect. 4.4 to estimate the mass accre- 
tion rate Mace (Fig- 8). The resolved double-peaked profiles of 
the allowed Fe i and Fe n lines in the spectrum of 08576nr292 
indicate their origin in a Keplerian rotating circumstellar disk. 
Furthermore, the CO-bandhead feature at 2.3 fim is likely also 
produced in the disk. This feature is a superposition of double 
peaked lines. Their peak separation is determined by the mass 
of the central object and the inclination angle, /, of the system, 
defined as the angle between the disk rotation axis and the line 
of sight. It can thus be used to estimate these parameters interde- 
pendently (see Sect. 6.1). For a detailed description of the spec- 
trum of 08576nr292 and a reconstruction of the system's geom- 
etry, we refer to EUerbroek et al. (201 1). 

The emission spectrum of 08576nr480 is dominated by lines 
from the cloud and the jet (see Sects. 3.1 and 3.2). The excep- 
tions are a few H i and O i lines and the Ca 11 infrared triplet, 
which are thought to originate in the accretion colunms and can 
be used to estimate the mass accretion rate; see Sect. 4.4. The 
CO-bandhead feature at 2.3 fim is also detected. An estimate for 
/ is not obtained from this feature due to insufl^icient signal-to- 
noise and spectral resolution. 

The position of the central source on the slit was determined 
with a gaussian fit on the spatial profile integrated over a spectral 
region surrounding the emission line. The continuum emission 
of 08576nr480 is not detected at /I < 1 .5 fim. For the lines in this 
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HH 1042 

[Olll]500 [01] 630 [Nil] 654 Ha 656 [SI I] 673 [Call] 729 [Oil] 731 [Mill] 737 01844 [Sill] 953 Hel 1083 [Fell] 1643 H2212I Bry2165 

! 



S 2 




-300 300 
v_sys (km s"') 



HH 1043 



[ONI] 500 [01] 630 [Nil] 654 Ha 656 [SI I] 673 [Call] 729 [Oil] 731 [Mill] 737 01844 [SHI] 953 Hel 1083 [Fell] 1643 H2212I Bry2165 




-300 300 
v_sys (km s"^) 



Figure!. Position- velocity diagrams of HH 1042 (top) and HH 1043 (bottom) for various lines, as labeled. The absolute flux scale is logarithmic. 
The underlying stellar continuum (at distance = 0''), where present, was subtracted using a gaussian fit. The ambient cloud produces emission in 
most lines at zero velocity (by definition) over the full length of the slit. In most lines the blue lobe of the HH 1042 jet is very prominent, while 
the red lobe suff'ers from extinction. The measured radial velocities in HH 1043 are significantly lower than those in HH 1042. 



HH 1042 

-200 -100 100 



[Fe II] 1643 nm 



HH 1043 




300 



100 

Vsvs (km s"') 



300 



Figure 3. Position- velocity diagrams of the [Fe ii] 1643 nm line of HH 1042 (left) and HH 1043 (right). The underlying stellar continuum (at 
distance = 0'') was subtracted using a gaussian fit. The positions of the knots are indicated. The dashed line indicates the position of the continuum 
source; km s~^ corresponds to the systemic velocity (see text). The remnant emission between -30 up to 50 km s~^ is a residual of the subtraction 
of a telluric OH emission line. 
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Figure 4. Sections of the continuum-normahzed on-source spectra of 08576nr292 (HH 1042, upper black hne) and 08576nr480 (HH 1043, lower 
gray line). Top: /-band; note the prominent contributions from the disk and accretion columns in 08576nr292, and the strong cloud emission in 
08576nr480 (e.g. the extensive H i Paschen series). Bottom: A'-band; both sources show prominent CO bandhead emission, likely produced by a 
Keplerian rotating disk. 



Table 2. Identified emission lines originating in the cloud and in the jet. 



A (nm)^ 


Ion 


Multiplet 


Cloud 


Jet 


A (nm) 


Ion 


Multiplet 


Cloud 


Jet 


A (nm) 


Ion 


Multiplet 


Cloud 


Jet 


372.603 


[On] 


IF 


+ 




738.818 


[Fen] 


14F 






+ 


1278.776 


[Fen] 


a^D-a^D 




+ 


372.882 


[On] 


IF 


+ 




745.254 


[Fe n] 


14F 






+ 


1281.808 


Hi 


5-3 


+ 


+ 


434.046 


Hi 


5-2 


+ 




763.754 


[Fe n] 


IF 






+ 


1294.269 


[Fe II] 


a^D-a^D 




+ 


486.133 


Hi 


4-2 


+ 




844.636 


Oi 


4 




+ 


+ 


1297.773 


[Fe II] 


a^D-a^D 




+ 


495.891 


[0 III] 


IF 


+ 




859.839 


Hi 


14-3 




+ 


+ 


1316.485 


Oi 


3P-3S° 


+ 


+ 


500.684 


[0 III] 


IF 


+ 




861.695 


[Fen] 


13F 






+ 


1320.554 


[Feii] 


a^D-a^D 




+ 


587.566 


He I 


11 


+ 




875.047 


Hi 


12-3 




+ 


+ 


1327.777 


[Fen] 


a^D-a^D 




+ 


630.030 


[0 I] 


IF 


+ 


+ 


886.278 


Hi 


11-3 




+ 


+ 


1533.471 


[Fen] 


a^F-a^D 




+ 


631.206 


[Sni] 


3F 


+ 




889.191 


[Fen] 


13F 






+ 


1599.472 


[Fen] 


a^F-a^D 




+ 


636.378 


[Oi] 


IF 


+ 


+ 


901.491 


Hi 


10-3 




+ 




1643.549 


[Fen] 


a^F-a^D 




+ 


654.804 


[Nn] 


IF 


+ 


+ 


905.195 


[Fen] 


13F 






+ 


1663.766 


[Fen] 


a^F-a^D 




+ 


656.280 


Hi 


3-2 


+ 


+ 


906.860 


[Sni] 


IF 




+ 


+ 


1676.876 


[Fen] 


a^F-a^D 




+ 


658.345 


[Nn] 


IF 


+ 


+ 


922.662 


[Fen] 


13F 






+ 


1680.652 


Hi 


11^ 


+ 


+ 


667.815 


He I 


46 


+ 




922.901 


Hi 


9-3 




+ 


+ 


1711.127 


[Fen] 


a^F-a^D 




+ 


671.644 


[S n] 


2F 


+ 


+ 


926.756 


[Fe n] 


13F 






+ 


1736.211 


Hi 


10^ 


+ 


+ 


673.082 


[S n] 


2F 


+ 


+ 


953.110 


[S ni] 


IF 




+ 


+ 


1744.935 


[Fe n] 


a^F-a^D 




+ 


706.525 


He I 


10 


+ 




954.597 


Hi 


8-3 




+ 


+ 


1797.103 


[Fe n] 


a^F-a^D 




+ 


713.579 


[Ar III] 


IF 


+ 




985.026 


[C I] 


IF 




+ 


w 


1809.394 


[Fe n] 


a^F-a^D 




+ 


715.516 


[Fe II] 


14F 




+ 


1004.937 


Hi 


7-3 




+ 


+ 


1817.412 


Hi 


9-4 


+ 


+ 


717.200 


[Fe II] 


14F 




+ 


1028.673 


[S II] 


3F 




+ 


+ 


1875.101 


Hi 


5-3 


+ 


+ 


725.445 


Oi 


20 


+ 




1032.049 


[S II] 


3F 




+ 


+ 


1895.310 


[Fe n] 


a^F-a^D 




+ 


728.135 


He I 


45 


+ 




1033.641 


[Sn] 


3F 




+ 


+ 


1938.770 


[Ni n] 


4p_2p 




+ 


729.147 


[Ca II] 


IF 




+ 


1037.049 


[Sn] 


3F 




+ 


+ 


1944.556 


Hi 


8-4 


+ 


+ 


731.992 


[On] 


2F 


+ 


+ 


1083.034 


He I 


1 




+ 


+ 


2058.130 


He I 




+ 


+ 


732.389 


[Ca II] 


IF 




+ 


1093.810 


Hi 


6-3 




+ 


+ 


2121.257 


H2 


1-0 8(1) 


+ 




732.967 


[On] 


2F 


+ 


+ 


1188.285 


[Pn] 




D2 




+ 


2165.529 


Hi 


7^ 


+ 


+ 


733.073 


[On] 


2F 


+ 


+ 


1256.680 


[Fen] 


a^D- 


a^D 




+ 


2222.685 


H2 


1-0 8(0) 
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Figures. Merged position- velocity diagrams of three individual lines traced with colors. Top: HH 1042; the stellar continuum of 08576nr292 
(at 0'') was subtracted using a gaussian fit. Bottom: HH 1043; no continuum removal was performed. Note the prominent emission of high 
excitation species (He i, [S iii]) in the shock regions where the velocity drops. The [O i] line peaks on- source, where the ejection mechanism 
operates, and in the shock regions. 
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Figure 6. The observed (i.e. not corrected for extinction) integrated flux of selected lines in the jet velocity range: -300 < v < -75 km s ^ (blue) 
and 75 < V < 300 km s"^ (red) for HH 1042; -300 < u < -30 km s"^ (blue) and 30 < u < 300 km s'^ (red) for HH 1043. Note the overlap of 
the two slit positions; the on-source emission profiles in the blue and red lobes trace different parts of the flow. Only emission above the average 
3cr noise level is shown. Note the different locations of the peaks and trends in intensity for e.g. the [Fe ii] and [He i] lines. 



spectral region we have used the source position derived from 
the ^-band continuum. In the plots of the jet spectra (Figs. 2, 
3, and 5) the stellar continuum (where present) was removed by 
subtracting the gaussian fit from the spatial profile. 

3.3. Emission from the jet 

In the spectra of both HH 1042 and HH 1043 we detect emis- 
sion at high blue- and redshifted velocities (100 - 220 km s"^) 
in opposite directions with respect to the source position, in a 
number of typical jet tracers (e.g. the [O i], [S ii] and [N ii] for- 
bidden lines in the optical; Fig. 2, Tab. 2). This emission traces 
the two lobes of bipolar jets that were already apparent from the 
SINFONI velocity map (Ellerbroek et al., 201 1). 

The emission line position- velocity diagrams (Figs. 2, 3 and 
5) show a velocity structure typical of jets, with successive ve- 
locity jumps of several tens of km s"^ as commonly observed 
in shocks. The observed lines are mainly from low ionization 
species, as is observed in jets from low-mass stars where the typ- 
ical shock velocities are - 30-40 km s"^ (e.g. Hartigan et al., 
1994, 2001). The terminal bow-shocks of jets can have much 
higher shock velocities, resulting in emission from high excita- 
tion species (e.g. [S iii]. He i). This is also seen in the shocks in 
our observations (e.g., knot E in HH 1042, Fig. 5). 

The three-color position- velocity diagrams (Fig. 5), pre- 
sented for the first time in the analysis of HH jets, highlight 
the variation of the excitation conditions along the jet. As usu- 
ally observed in jets, strong emission is present in transitions 



of refractory ion species such as Fe, Ca and Ni, which are de- 
pleted onto dust grains in the ISM. When dust evaporates in the 
jet launching region and/or in shocks, these ion species are sub- 
sequently released in the gas phase. In the following, the main 
morphologic and kinematic characteristics are discussed for the 
two sources individually. 



3.3.1. HH 1042 

The emission from the bipolar jet HH 1042 in the bright 
[Fe ii] 1643 nm line is shown in Figs. 1 and 3. The blue lobe 
includes seven knots labeled A to G, and extends up to 13'' from 
the central source. It is covered by the slit up to knot F, at ~ 9'' 
from the source, then it terminates in a non-collimated structure, 
knot G, which is visible in the SINFONI map (Fig. 1). The red 
lobe extends up to 9'' from the source. Beyond this point it is not 
detectable, most likely due to foreground extinction. The five 
knots in the red lobe are labeled A' to E' (Figs. 1 and 3). The jet 
produces emission lines in allowed transitions of H, He and O, 
and forbidden transitions of O, P, S, Ca, Fe and Ni (Tab. 2). 

Fig. 2 shows the emission along the jet in a selection of lines. 
We see that the strongest emission is detected in the [Fe ii] lines, 
while also H i. He i and [S iii] are prominent along the jet. The 
maximum brightness is at the position of the bright knot E with a 
flux in the [Feii] 1643 nm line of 2.05 ±0.01 xlO'^'^ ergs"^ cm'^. 
The line fluxes measured in each knot are listed in Tables A.l 
and A.2 in the Appendix. The root-mean- square errors were cal- 
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culated from the error spectrum (based on readout noise, flat- 
field, dark and bias) provided by the X-shooter pipeline. 

The composite position-velocity diagrams in Fig. 5 show 
variations of the excitation conditions along the jet. [O i] emis- 
sion dominates on- source, while moving along the jet an in- 
creasing degree of ionization is seen. High ionization/excitation 
lines, e.g.. He i and [S iii], dominate in the tail end, particularly 
knots E and F. This could be due to the terminal shocks being 
the strongest in the jet. 

The line flux along the jet (Fig. 6) increases significantly 
beyond knot E in the blue lobe; in the red lobe, the emission 
suff'ers from extinction, which increases dramatically beyond 5'' 
from the source (see Sect. 4.1). The velocity in the jet is approx- 
imately 130 km s"^ at the base, then increases up to 220 km s"^ 
right before knot E, after which it falls back to 140 km s"^ in the 
blue lobe (see Fig. 9). Similarly in the red lobe, the velocities 
vary between 120 km s"^ and 210 km s"^ . 

Note that no clear correlation exists between the knots in 
the blue (A-G) and the red (A'-E') lobes in terms of posi- 
tion and velocity. Asynmietries in velocity between the blue 
and red lobe are commonly observed in jets; this is usually at- 
tributed to an interaction with the ISM (see e.g. Hirth et al., 1994; 
Melnikov et al., 2009; Podio et al., 2011). However, in the case 
of HH 1042, the average velocity is roughly equal in both lobes, 
although there is an uncertainty in the value of fsys- The variation 
in line flux and velocity along the jet on either side of the source 
is somewhat symmetric, although a definitive match between the 
knots in the blue and the red lobes cannot be made. We further 
conmient on this in Sect. 5, where the (a)symmetry of the jet is 
compared with models for the outflow. 

3.3.2. HH 1043 

The emission knots of HH 1043 are more clearly separated than 
those in HH 1042 (Figs. 1 and 3). The blue lobe consists of two 
knots (A and B), the latter of which ends in a bow-shock shape 
that is spatially resolved on the SINFONI linemap (Fig. 1). The 
red lobe is separated into four knots named from A' to D'. The 
blue and red lobes both have projected lengths of 8.5". Note that, 
as in HH 1042, there is no clear symmetry between the blue and 
red lobes in terms of the positions of the knots. The same lines 
that are detected in HH 1042 are also present in HH 1043, with 
the addition of several H2 lines and higher H i transitions (up to 
the Balmer, Paschen and Brackett jumps) due to cloud emission. 

The brightest lines are again the [Fe 11], Ha, He i and [S ni] 
lines (Figs. 2 and 6). The terminal knot B has an integrated flux 
of 1.03 ± 0.01 X 10-1^ erg s"^ cm'^ for the [Fe 11] 1643 nm line. 
The line fluxes measured in each knot are listed in Table A. 3 in 
the Appendix. Like in the other HH object, the on- source knot 
A has the strongest [O i] flux, while in the terminal bow-shock 
(knot B), [O i], [S II], [S III] and He i are bright (Fig. 5). 

The resolved bow-shock shaped feature at the end of the blue 
lobe and the comparable brightness of the two lobes suggest that 
the inclination of the jet is quite high. The average velocities in 
both lobes are somewhat asymmetric. Some knots in the difl'er- 
ent lobes can be matched: knots B and D' are located at ~ 8'' 
either side of the source, while knots A and B' are at diff'erent 
distances (T' and respectively). Knot A' has no visible coun- 
terpart in the blue lobe. 

4. Physical properties of the jets 

The detected emission lines contain information on the physical 
conditions of the gas in the jet. By using selected line ratios as 



diagnostic tools one can estimate the electron density and total 
density, the ionization fraction and the temperature. In particu- 
lar, the observed line ratios can be compared: (i) with those pre- 
dicted by shock models (e.g. Hartigan et al., 1994); or (ii) with 
ratios computed assuming that the employed forbidden lines are 
optically thin and collisionally excited, i.e. assuming that the in- 
teraction with the radiation field is negligible (e.g. Bacciotti & 
Eisloff'el, 1999; Podio et al., 2006). In Table 3 we summarize the 
line ratios used in this paper to estimate the physical properties 
of HH 1042 and HH 1043. Furthermore, in Sect. 6, Fig. 13, some 
ratios are compared to those observed in similar sources. 

In order to improve the signal-to-noise ratio, the fluxes are 
integrated spectrally over their profile and spatially over the de- 
fined knots. Tables A.l, A.2 and A.3 display the integrated flux 
for all the emission lines in the knots where the flux exceeds the 
background noise by a factor 3. Table 3 lists the values of the 
physical properties derived from selected line ratios. 

4.1. Extinction 

The [Fe 11] 1643/1256 nm and 1643/1321 nm line ratios only de- 
pend on the intrinsic ratio of Einstein coefficients because the 
considered lines share the same upper level. Thus, the difl'erence 
between observed and theoretical [Fe 11] ratios is a direct tracer 
of the visual extinction (Ay). However, the values for Ay inferred 
from the [Fe 11] 1643/1321 nm line ratio are systematically lower 
than those estimated from the [Fe 11] 1643/1256 nm line ratio, 
because of the uncertainties afl'ecting the computed Einstein co- 
efficients. Moreover, when using different sets of Einstein coeffi- 
cients in the literature (e.g. Nussbaumer & Storey, 1988; Quinet 
et al., 1996), diff'erent values for Ay are obtained. 

This issue is discussed in some detail in Nisini et al. (2005), 
Podio et al. (2006) and Giannini et al. (2008). In particular, 
the latter authors compared the Ay values obtained from line 
ratios assuming diff'erent sets of Einstein coefficients with Ay 
values estimated with other, independent methods and showed 
that the most reliable estimate of Ay is obtained by using the 
[Fe n] 1643/1256 nm line ratio and the Einstein coefficients by 
Quinet et al. (1996). We have adopted the same approach in this 
work. The [Fe 11] 1643/1256 nm line ratios and the obtained Ay 
values are shown in Fig. 7 and summarized in Table 3. A global 
uncertainty on these values may be caused by the value of the 
total- to- selective extinction Ry. Throughout this paper we adopt 
the average Galactic value of Ry = 3.1 and the extinction law 
from Cardelli etal. (1989). 

In HH 1042, the red lobe appears much fainter than the blue 
lobe, which can be explained by the extinction trend in Fig. 7. 
This is consistent with the red lobe disappearing into (or behind) 
the molecular cloud, as was proposed in EUerbroek et al. (201 1). 
Even when corrected for extinction, the flux level in the red lobe 
is 2 - 3 times less than in the blue lobe. 

The on- source value of Ay measured from the [Fe 11] 
1643/1256 nm ratio is much smaller (Ay = 0.79 ± 0.21) than 
that estimated by EUerbroek et al. (2011) from fitting the SED 
to a disk slope (Ay = 8 ± 1). The measured on-source line ratio 
might be closer to unity than its true value as a result of resid- 
uals of the continuum subtraction, leading to an underestimate 
of the true extinction. However, even within this uncertainty, the 
on-source extinction would still be much lower than the value 
derived from the SED-fitting. This suggests that between the jet, 
traced by the [Fe 11] lines, and the protostar a dusty shell or the 
disk might further obscure the photosphere. This phenomenon is 
not uncommon; for example, the extinction towards DG Tau B 
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Figure?. Top panels: [Fe ii] 1643 nm and 1256 nm line fluxes as a function of position. Bottom panels: [Fe ii] 1643/1256 nm flux ratio, integrated 
over knots. The right hand axis indicates the corresponding values of Ay (Cardelli et al., 1989; Quinet et al., 1996). Left: HH 1042; Right: HH 1043. 



is much higher than the extinction derived from its jet emission 
lines (Podio et al, 2011). 

In HH 1043, somewhat higher extinction values are mea- 
sured in the red lobe than in the blue lobe, consistent with the 
slightly inclined position of the jet as discussed in Sect. 3.3.2. 

4.2. Electron density, temperature and ionization fraction 

The values for the jet physical conditions derived from the ob- 
served line ratios are listed in Table 3. The electron density 
Hq is estimated using the [S ii] 673/671 nm and the [Fe ii] 
1643/1533 nm line ratios. These diagnostics yield values that 
agree to within a factor two for most knots in the blue lobes. 
However, in the red lobes of both objects, the values de- 
rived from [Fe ii] lines are significantly higher than those esti- 
mated from [S ii] lines. This may be because they trace a zone of 
the post-shock cooling region which is located further from the 
shock front where the gas is more compressed, as discussed by 
Nisini et al. (2005) and Podio et al. (2006). 

The electron temperature Tq is calculated from the [Fe ii] 
1643/861 nm ratio, which is independent of the electron density 
to within an order of magnitude of (Nisini et al., 2005). We 
detect a decreasing trend in electron temperature moving away 
from the sources. This is commonly observed in HH jets; see e.g. 
Bacciotti & Eisloffel (1999); Podio et al. (2006, 2011). 

Predictions by shock models (Hartigan et al., 1994), as well 
as theoretical line ratios demonstrate that, once Hq is deter- 
mined, the [N n]/[0 i] line ratio is almost independent from Tq, 
hence can be used to estimate the hydrogen ionization fraction 
Xq = nu+/nu- It should be noted that the theoretical line ratios 
are estimated by assuming that lines are collisionally excited and 
that charge exchange between O, N and H is the dominant pro- 
cess determining the hydrogen ionization fraction. 

We note that [O i] emission in HH 1042 is only detected on- 
source and in the bright knot E, while [N ii] is detected all along 
the jet, but very weakly on-source (see Figs. 2 and 5). Thus, only 
in knot E can we can compute the [N ii]/[0 1] line ratio and derive 
an estimate of Xq from the Hartigan et al. (1994) shock models. 
We indeed find a high ionization fraction in knot E, where the 
steepest velocity gradient is located, indicating that shocks may 
contribute to the increased ionization conditions. From the upper 
limit on the on-source [N ii]/[0 1] line ratios we derive Xq < 0.025 
in knot A. In HH 1043, the [O i] and [N ii] emission is too weak 
(less than 3cr) along the whole jet, making a reliable estimate of 
Xq impossible. 



The [O i] emission peaks on-source, where it is formed in the 
energetic (disk) wind which constitutes the base of the jet (Cabrit 
et al., 1990). The [O i] velocity in knot A of HH 1042 coincides 
with the blueshifted absorption component of the Ca ii triplet 
lines (see Fig. 4), strongly suggesting that these lines originate 
in the same medium. Finally, by using our Hq and Xq estimates, 
and assuming that the free electrons are due to the ionization 
of hydrogen atoms, we derive an estimate of the total density 
nn = Hq/Xq = 6.46!|-^^ x 10^ cm"^ in knot E of HH 1042. 

The uncertainties on these estimated physical quantities are 
dominated by diff'erent eff'ects. The uncertainty on Hq as derived 
from the [S ii] and [Fe ii] ratios is dominated by the error on the 
line fluxes, as both pairs of wavelengths are close together, mak- 
ing the efl'ect of extinction negligible. The uncertainty on Tq is 
dominated by the error on Ay as the [Fe ii] 1643 nm and 861 nm 
lines used for that estimate are further apart in the spectrum. The 
uncertainty in Xq is dependent on the errors in Ay and the line 
flux in equal measure. 

4.3. Mass outflow rate 

An important quantity in jet dynamics is the mass outflow rate, 
Mjet. It determines how much mass and linear momentum is 
injected in the surrounding cloud, and when the jet rotation 
is known, how much angular momentum is removed from the 
YSO. The ratio of the mass outflow rate to the mass accretion 
rate. Mace, determines the eflftciency of the star formation pro- 
cess. Magneto-hydrodynamic models of jets typically adopt val- 
ues in the range MjQt/M^cc ~ 0.01-0.1 (Konigl & Pudritz, 2000; 
Shu et al, 2000; Cabrit, 2009). 

The mass outflow rate can be estimated from the observed 
line fluxes and their ratios using three difl'erent methods (see 
Tab. 3): 

(i) By multiplying the total density nu with the transverse cross- 
section 7r7?j and the deprojected velocity | f;j | = v_i/ cos / of 
the jet (Podio et al., 2006): 

Mjet = /dmnnuTTR^lvjl (1) 

where mn is the proton mass and ji = 1.24 was adopted 
for the mean atomic weight. For both objects, we adopted 
a value of R] = 200 AU for the jet radius. This is the av- 
erage measured half width at half maximum of the spatially 
resolved [S ii] intensity profile in similar HH objects, which 
ranges from 75 - 300 AU (e.g. Mundt et al, 1991; Reipurth 
et al., 2000, 2002; WasseU et al, 2006). The inclination is 
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Table 3. Physical parameters and mass loss rates estimated in the brightest knots along the HH 1042 and HH 1043 jets. 



HH 1042 HH 1043 

Blue lobe Red lobe Blue lobe Red lobe 



Quantity 


Diagnostic (nm) 


Kei. 


knot A 


knot E 


knot 


knot A 


knot B 


knot A' 


knot 


Ay (mag) 


[Fe ii] 1643/1256 


Q96 


0.79!0-2| 


•^•^"-0^09 


-'"^ -0.34 


c o /I +0 20 
^••^^-0.20 


A ci +0 14 


c /I +0 22 
-^•^^-0.22 


/c niz+0 37 

6.76:0-37 


Hq (ID cm ) 


\p IlJ o/J/o/i 


(JrUo 


> o./U 




1 S<d+0-32 
-"^•-^ -0.27 


A ao+2.17 


on+1.23 

6.90_o97 


c /I7+1.88 






[Fe n] 1643/1533 


N05 


9 98+3-^'^ 

^•^°-0.79 


3 51+0-23 

■^-^^-0.23 


7 1 O+0.58 
/•^^-0.60 


8.04!0-45 


Q 07+0.40 

^•^/-0.41 


28.6!|o 


20.3!^;^ 


Te (103 K) 


[Fen] 1643/861 


N05 


10 6+0-3 


5 07+0-03 


4 7g+o.ii 


6 07+0-08 

'^•^ -0.08 


c 9Q+0.05 
^•^^-0.05 


c 77+0.07 
^•^ ^-0.07 




Xq 


[N n] 654 /[O i] 630 


H94 


< 0.025 


-0.7 












Mje, (10-8 yj.-l) 




P06 




Q ^Q+2.02 
^•-^^-2.08 














^[S II] . 


H95 




(> 0.09) 


(> 0.07) 


(> 0.05) 


(> 0.11) 


(> 0.10) 






^[O I], 


H95 




(> 0.03) 


(> 0.05) 


(> 0.02) 


(> 0.02) 


(> 0.02) 






^[O I]» T'e, I'shock 


KT95 


5.68!3-| 


1 o Q+10.7 


10.3:5-9 


09 1+18.4 
■^^•^-18.4 


27 1+^5-5 







H94: Hartigan et al. (1994); H95: Hartigan et al. (1995); KT95: Kwan & Tademaru (1995); N05: Nisini et al. (2005); OF06: Osterbrock & Ferland (2006); P06: Podio et al. (2006); 
Q96: Quinet et al. (1996) 



17.8°_+io4^ for HH 1042 and i 



60°_\l^ for 



(ii) 



estimated at / 
HH 1043 (see Sect. 6.1) 
From the [O i] and [S ii] emission line luminosities, if it 
can be assumed that all oxygen is neutral and sulphur is all 
singly ionized. The total line luminosity is then proportional 
to the number of emitting atoms in the observed volume. 
Adopting coUisional coefficients and critical densities ricr as 
in Hartigan et al. (1995, equations A8-A10), we have: 



jet 



and 



9.61 X 10"^ 1 + 



1 + 



^cr,[OI] 



Mjet = 1.43 X IQ-^ 



^cr,[SII] 

rie 



HOI] 



^knot 



^knot 



Moyr-i (2) 



Moyr 



(3) 

(iii) From the [O i] line luminosity if this is produced by post- 
shock cooling, by using Eq. (A 14) from Hartigan et al. 
(1995), which is adopted from Kwan & Tademaru (1995): 



^j^^ " — ■ 



(4) 



We are unable to determine the shock velocity f;sh from spa- 
tially resolved shock fronts. Typical values from the models 
discussed in Hartigan et al. (1994) are in the range i^jet/f^sh ~ 
5 - 10; we adopt f^jet/f^sh = 10- It is assumed that a fraction 
1 // of the total luminosity radiated below 6000 K is in the 
[O i] 630 nm line. As in Kwan & Tademaru (1995), we adopt 
/ = 3.5. 

The uncertainties on the Mjet estimates from method (i) are 
mainly due to the assumption of the jet radius and the uncertainty 
on the inclination angle, while for method (ii) and (iii), which 
are based on the luminosities of forbidden lines, the main source 
of error is the uncertainty on the extinction and the distance to 
the source and - for method (iii) - the estimated shock velocity. 
An uncertainty in the absolute flux calibration due to the cor- 
rection for slit losses may be present, although this marginally 
contributes to the error budget. The excitation models used in 
method (ii) assume that all oxygen is neutral and all sulphur is 
singly ionized. As emission in [O n] and [S iii] lines is detected 
in the jets, the estimates for Mjet obtained from this method are 



considered lower limits. The high excitation conditions in the 
jet may point to an external source of radiation (see Sect. 6.2). 
However, since even higher ionization species are detected in the 
cloud and not in the jet, any external radiation field causing this 
emission is not expected to aff'ect the estimates on the physical 
conditions and mass outflow rates in a significant way. 

Method (i) can only be applied in knot E of HH 1042, be- 
cause it is the only knot where we can retrieve an estimate of 
Xq. Method (iii) can be applied in those knots where [O i] is de- 
tected and we have an estimate of Ay. On-source, where the gas 
is thought to be almost neutral as indicated by the derived up- 
per limit on for HH 1042 (see Tab. 3), we cannot apply this 
method because we do not have a reliable estimate of n^, due to 
the faintness of the [S n] lines. 

Tab. 3 summarizes the values and lower limits obtained by 
applying the explained methods. The absolute values of Mjet ob- 
tained by methods (i) and (iii) agree well and indicate Mjet ~ 
10-^ Mo yr-i in both HH 1042 and HH 1043 (see Tabs. 3 and 5). 
In both jets, no significant asymmetry is found between the red 
and blue lobe values of Mjet. 
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Figure 8. Values for Mace of 08576nr292 derived with Ay = 0.79 + 0.21 
(open circles) and Ay = 8 ± 1 (closed circles). 
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Table 4. Estimated mass accretion rate. 



Line, A 


log Flux ' 


' log Flux 


log Mace 


Reference^ 


(nm) 


(erg cm 


(erg cm~^) 


(Mo yr-i) 








08576nr292 (HH 1042) 








Av = 


Av = S + l 


Ay = 8 + 1 




Ca II K 393 


-15.28 


-10.66 + 0.04 


-5.95 ± 0.30 


HH08 


H7434 


-15.71 


-11.43 ±0.06 


-6.03 ± 0.21 


HH08 


1^486 


-14.83 


-11.05 ±0.04 


-5.95 ±0.21 


HH08 


He 1 501 


-15.21 


-11.59 ±0.06 


-5.66 ±0.31 


HH08 


He 1 587 


-15.35 


-12.46 ±0.17 


-5.51 ±0.74 


HH08 


NaiD589 


-14.96 


-12.03 ±0.07 


-6.20 ± 0.73 


HH08 


[0 i] 630 


-14.36 


-11.64 ±0.04 


-6.09 ± 0.84 


HH08 


Ha 656 


-13.02 


-10.42 ±0.01 


-5.77 ± 0.40 


HH08 


Ca II 854 


-12.89 


-11.19 ±0.01 


-6.02 ±0.51 


HH08 


Pay 1093 


-13.14 


-12.02 + 0.02 


-5.89 + 0.62 


G08 


PaJS 1281 


-12.68 


-11.82 + 0.01 


-6.20 + 0.62 


M98 


Bry 2165 


-12.87 


-12.50 + 0.01 


-5.97 + 0.86 


M98 


Ua 10% width 


686 km s" 


-1 


-6.24 + 0.30 


N04 






08576nr480 (HH 1043) 








Av = 


Ay = 12 ±3 


Ay = 12 + 3 




Can 854 


-15.66 


-13.11 ±0.09 


-7.97 ± 0.60 


HH08 


Par 1093 


-14.32 


-12.64 ±0.05 


-6.73 ± 0.65 


G08 


PayS 1281 


-13.82 


-12.52 ±0.03 


-7.00 ± 0.67 


M98 


Bry 2165 


-13.96 


-13.40 + 0.03 


-7.11+0.93 


M98 



The errors are less than 0.01 dex. 

G08: Gatti et al. (2008); HH08: Herczeg & Hillenbrand (2008); M98: Muzerolle et al. 



(1998a); N04: Natta et al. (2004) 



Table 5. Mass outflow and accretion rate. 



Object 


HH 1042 / 
08576nr292 


HH 1043 / 
08576nr480 


Mjet (Mo yr-1) 
M^cc (Mo yr^) 


8.86:j g X 10"^ 

i.io:o-29 X 10-6 


2.92:};!^ X 10"^ 
5.50!|g X 10-^ 



4.4. Accretion rate 

The accretion rate, Mace, is derived from the accretion luminos- 
ity, Lace, which can be determined by measuring the UV excess 
flux emitted by the accretion flow close to the star (e.g. Hartigan 
et al., 1995; Gullbring et al., 1998). Radiative transfer models 
and spectra show that certain spectral lines are formed in the 
accretion flow (Hartmann et al., 1994). Therefore one would ex- 
pect that the line strength correlates with the accretion luminos- 
ity, which was confirmed for CTTS in subsequent studies (e.g. 
Muzerolle et al., 1998a,b). This correlation between line strength 
and Lace is consistent across the mass spectrum, from brown 
dwarfs up to HAe stars (e.g., Natta et al. 2004; Mendigutia et al. 
2011, see also Fig. 14). 

The mass accretion rate can be related to the accretion lumi- 
nosity, following Gullbring et al. (1998). The accretion luminos- 
ity is equal to the amount of energy per unit time released from 
the gravitational field when material falls onto the stellar surface 
- along magnetic field lines - from the radius Rin where the disk 
is truncated by the stellar magnetic field: 



I ^ ^* \ Lacc^* 



(5) 



We adopt M^/R^ ~ 1 (in solar units) for the central star, consis- 
tent with pre-main sequence models (PMS, Siess et al., 2000), 
and assume a typical value of Ri^ ~ 5 (Shu et al, 1994). 

Fig. 8 shows that for 08576nr292, adopting an extinction 
Ay = 0.79 ± 0.21 found from the [Fe 11] line ratio, the accre- 
tion diagnostics are not consistent. Instead, using Ay = 8.0+1.0 



as estimated from SED fitting by EUerbroek et al. (201 1), we ob- 
tain consistent results for tracers across the entire spectral range: 
log Lacc/^o = 1 .53 ± 0. 10, which is of the same order as the stel- 
lar luminosity (Sect. 6). Subsequently log Mace = -5.96 ± 0.10 
Mo yr"^. This confirms that the on-source extinction is much 
higher than that estimated at the base of the jet from [Fe 11] lines. 

The spectrum of the driving source of HH 1043, 08576nr480, 
only exhibits a few accretion tracers (Ca 11, Bry, Pa^ and Pay). 
Using these tracers, the accretion luminosity is found to be 
logLacc/Lo = 0.23 ± 0.34, and log Mace = -7.26 ± 0.34 M© 
yr~^. We adopted Ay = 12 ± 3 from the SED-fitting method de- 
scribed in EUerbroek et al. (2011). Tab. 5 lists the derived mass 
outflow and accretion rates. 

5. Kinematics analysis 

In this section, we present our analysis of the kinematics of the 
HH 1042 jet, as well as an interpretative model that we use 
to simulate the recent outflow history of the jet. Our aim is to 
learn more about the central YSO and the jet launching mech- 
anism from the fossil record of the outflow history contained 
in the jet. We focus on the position- velocity diagram of the 
[Fe 11] 1643 nm line. This line is one of the brightest in the spec- 
tra and traces the largest velocity range. An additional advantage 
is that the detector resolution is the highest in the near-infrared 
arm (Av ~ 26 km s~\ see Table 1). Note, however, that given the 
spatial and spectral resolution of the observations, and the lack 
of proper motion measurements, one cannot formulate a unique 
model that reproduces the data. Instead, in this section we simu- 
late the general shape of the emission pattern, derive timescales 
relevant for the ejection mechanism, and draw qualitative con- 
clusions about the physics within the flow. 

The position-velocity diagram of the [Fe 11] line shows an 
outflow variable in velocity (Figs. 3 and 9). This kinematic struc- 
ture is assumed to be the result of an outflow which varies in both 
velocity and mass outflow rate at the launch site at the base of 
the jet. The outflow rate is reconstructed by comparing the data 
to a ballistic outflow model that assumes a launch mechanism 
that is either stochastic or periodic. The simulation consists of 
two ingredients: a characterization of the temporal variation of 
the mass outflow rate at the base of the jet (the 'input physics') 
and a description of the flow of material through the jet (the 'in- 
teraction physics'). The former is reconstructed in Sect. 5.1 for 
both the stochastic and the periodic mechanisms; the latter is ex- 
plained in Sect. 5.2. The results of the simulations are presented 
in Sect. 5.3. The kinematics of HH 1043 are not simulated, as 
there are too few emission knots along the jet to constrain the 
model parameters. 

5.1. 'Input physics': outflow rates v(t), m(t) 

Time variability is a known property of accretion-ejection mech- 
anisms and has been measured in both the accretion luminosity 
(e.g. Herbst et al., 1994; Alencar & Batalha, 2002; Hillenbrand 
et al., 2012) and outflow activity (e.g. Micono et al., 1998) of 
YSOs on timescales ranging from days to years. It has been ar- 
gued by Hartmann & Kenyon (1985) that this reflects an intrin- 
sic variability of the accretion process. One may therefore ex- 
pect all observables correlated with accretion, i.e. the luminosity, 
mass accretion and outflow rates as well as outflow velocities, to 
be variable in time. The characteristics of these variations are 
not well constrained, but they are expected to be either purely 
periodic (e.g. due to disk rotation or binary interaction), quasi- 
periodic (e.g. due to the interplay between magnetic stress and 
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Figure 9. Top: Peak radial velocity, |i;|cos/, and FWHM of the velocity components detected along the HH 1042 jet, obtained by means of a 
multiple gaussian fit of the [Fe 11] 1643 nm line velocity profile at each spatial pixel (0.2'') along the jet. The size of the symbols represents the 
integrated flux of each velocity component. Bottom: Position- velocity diagram of the [Fe 11] 1643 nm line in HH 1042. The j-axis corresponds to 
the absolute value of the radial velocity, in order to show the (a) symmetry between the two lobes. The flux is corrected for extinction and the red 
lobe (left) is enhanced by a factor of 3, in order to compare the emission patterns of both lobes. The source is located at 0'' in position space. The 
letters indicate the knots which were integrated in flux for the kinematics analysis. The size of the circles corresponds to the integrated flux of the 
knots. 



pressure in the accretion disk), or stochastic (e.g. due to chaotic 
processes that depend on many physical parameters). Different 
periodicities and timescales, tracing different mechanisms, may 
exist within one accretion system. 

In principle, both the outflow velocity, v(t), and the mass 
outflow rate, m(t) can be variable. A varying v(t) results in dif- 
ferences in velocity and line flux (because of the formation of 
shocks) along the jet; a varying m(t) introduces a variation in 
density and hence line flux along the jet. In this section, we ex- 
plore how well we can reconstruct the outflow velocity and mass 
outflow rate from the data. 

We estimate the launch time of the material along the jet 
from its present position and velocity. From the position- velocity 
diagram of the [Fe 11] 1643 nm line, a one-dimensional spec- 
tral profile is extracted at every pixel of width = 0.2'' along 
the jet. To this emission profile, a superposition of one or more 
one-dimensional gaussian functions is fitted. We have made use 



of the IRAF routine splot, which deblends multiple gaussian 
components and calculates errors on the fit parameters based on 
a Poisson noise model. The number of components is increased 
until adding another component does not significantly improve 
the fit; in most cases this amounts to one or two components per 
row. For every spatial pixel row of every emission line, this re- 
sults in a list of emission components, their velocities, widths, 
and fluxes (Fig. 9, top). 

To increase the signal-to-noise ratio, these emission compo- 
nents are grouped into knots in position-velocity space (Fig. 9, 
bottom). The knots are similar to those defined in Sect. 3; some 
are split into several components with diflTerent radial velocities 
(labeled with subscripts, e.g. Ei, E2). The velocities and posi- 
tions of the knots are determined by averaging them over the 
values of the constituent components, where the velocities are 
weighted by their inverse squared errors. 
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Figure 10. Radial velocity, | v \ cos /, versus launch time, ^launch tan /, as- 
suming no collisions along the flow. Symbols represent the knots in the 
blue lobe (above, blue circles) and red lobe (below, red circles) as de- 
fined in Fig. 9 (bottom). The area of a circle corresponds to the line flux 
integrated over the knot. 



We estimate for every knot with measured line-of- sight ve- 
locity, f;cos/, and projected position in the plane of the sky, 
X sin /, the launch time, by assuming that no collisions occur in 
the flow (defined as the case f) = -I, see next section): 



^launch tan / 



xsmi 



vcos I 



(6) 



The results are displayed in Fig. 10. The intervals between the 
knots and the velocity variation appear to be somewhat regular, 
particularly in the blue lobe. One might suspect that this reflects 
a purely periodic outflow rate, with single-mode sine waves for 
v(t) and m(0, which would be the simplest case of an outflow 
rate which is variable in time. However, apparent regularities in 
time series may wrongly suggest periodicity to the human eye. 
A stochastic or quasi-periodic outflow rate may also produce 
an apparently periodic signal in the outflow pattern. We there- 
fore consider two possibilities for the variations in v(t) and m(t)\ 
a stochastic variation, the parameters of which are obtained by 
Fourier analysis; and a purely periodic variation, the parameters 
of which are obtained by an iterative, direct fitting procedure. In 
both cases, a diflTerent periodicity is allowed for v(t) and m(t). 
From this point on we refer to these two cases of simulated out- 
flow rates as 'stochastic model' and 'periodic model'. 

5.1.1. Stochastic model 

In order to see whether a periodic signal is present in the time 
series data displayed in Fig. 1 1 (top), we perform a Fourier anal- 
ysis to the (v cos /, ^launch tan /) coordinates and the (F, ^launch tan /) 
coordinates, where F is the flux per gaussian component. Only 
the blue lobe data are fitted, as the signal is highest there. 

Fig. 11 (bottom) shows the Lomb-Scargle periodogram 
(Scargle, 1982) for the velocity and flux time series. There is 
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Figure 11. Top: Timeseries (v, tiaunch) and (F, tiaunch) obtained by fitting 
the diff'erent Gaussian components in the emission profile (Fig. 9) as 
described in Sect. 5.1 and applying Eq. (6). Bottom: Fourier transform 
(Lomb-Scargle periodogram) of the (v, tiaunch) and (F, tiaunch) coordinates 
in the blue lobe. The best-fit power-laws (dashed lines), minus the white 
noise constant y, are overplotted. The dotted line indicates the largest 
period Ti tan / in the periodic model (Tab. 6); the flux periodogram (and 
the velocity periodogram, albeit very weakly) also peaks at this value. 



no clear indication for a periodic signal in the velocity data, 
but increasing power at longer periods indicate the presence of 
a noise component. The flux data shows a peak at cojidini ~ 
0.08 rad yr"^ (T tan / ~ 80 yr) suggesting a (quasi- )periodic pro- 
cess on long timescales. Incidentally, this period coincides with 
the timescale found for the periodic model (see Sect. 5.1.2). We 
fit both periodograms, using a Maximum Likelihood approach 
appropriate for power spectra (see e.g. Vaughan, 2005), with a 
power law plus a constant of the type 



P(v) = py-^ + y , 



(7) 



where a is the power law index, is a normalization term and y 
is a constant to account for the presence of white noise in the pe- 
riodogram. The velocity periodogram is well fit with this model, 
with a power-law index of = 1 .50+0.04. The flux periodogram 
is fitted by a power-law with index a = 1.98 ± 0.05 but for the 
peak at ~ 80 yr. 
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Table 6. Input parameters for periodic outflow rate, HH 1042. 



Parameter 


Value 


Parameter Value 


fo cos / ^ 


170kms-i 


01 0.27 (blue) 


Vi COS / 


39.8kms-i 


-0.38 (red)^ 


Ti tan/ 


83.0 yr 


02 0.18 (blue) 


72 tan/ 


12.0 yr 


-0.47 (red) 



The mean velocity vq cos i is fixed. 

The phase offset between the red and blue lobes is fixed at A0 = 0.65 in all simulations. 



In order to reconstruct the modulation in the outflow rate ac- 
cording to a noise process with the properties of periodograms 
for V and F, we employed the method by Timmer & Koenig 
(1995) to simulate time series from periodogram data. We sim- 
ulated 1000 time series v(t) and m(t) from the fit to the velocity 
and flux periodograms, with a time resolution in the simulated 
time series of 10 years (which roughly corresponds to the spatial 
resolution of the detector). These time series served as input for 
the model described in Sect. 5.2. 

5.1 .2. Periodic model 

An alternative approach to reconstructing the modulation of the 
outflow rates v(t), m(t) is by assuming them to be analytic, purely 
periodic functions: 



v(t) = Vq -\-v\ sin 27r 



m(t) = rho 



1 -Fsin27r| — 02 

^2 



(8) 
(9) 



Here, T\ and T2 are the periods of the oscillations in the velocity 
and mass outflow rate, respectively; 0i and 02 are the relative 
phases of the oscillations in the velocity and mass outflow rate, 
respectively. The parameters vq and vi are the mean and ampli- 
tude of the velocity variation; the parameter rho is scalable and is 
normalized to fit the flux level in the data. As no clear symmetry 
exists between the knots in the red and blue lobe, a phase off'set 
A0 was introduced between the blue and red lobe outflow veloc- 
ity. Summarizing, this amounts to seven free parameters for the 
periodic outflow rates: vq, vi, T\, T2, 0i, 02, and A0. 

A set of optimal values for these input parameters was ob- 
tained by fitting sinusoids to the data displayed in Fig. 10 to the 
knot velocities (resulting in v(t)) and to the injection intervals 
(resulting in m(t)). In order to fit the data, a phase off'set between 
the red and blue lobes was introduced by setting A0 = 0.65. The 
model parameters listed in Tab. 6 thus obtained, result in a sim- 
ulation that best represents the observed emission in position- 
velocity space. In the next section, a ballistic model is described 
with which we simulate the flow, with Eqs. (8) and (9) as the 
input velocity and mass outflow rate. 

5.2. 'Interaction physics': model for a ballistic flow 

In this section we describe a one-dimensional model for the en- 
ergy loss along a ballistic flow in a jet. This model is inspired 
by the approach of Raga et al. (1990), who solved the inviscid 
Burgers equation for a variable flow in one dimension. In later 
work (Raga et al., 2012) it was shown that the shock fronts (i.e. 
the knots) along such a flow merge by inelastic collisions. In our 
approach, we describe the flow in a Lagrangian picture, in terms 
of the collisions of discrete parcels of gas. As explained below. 



we parametrize between the limiting cases of a flow with inelas- 
tic collisions, and a free flow without collisions. 

In our simulations, a sequence of discrete parcels of varying 
mass jfii = mAt is ejected with a constant time interval At. The 
time interval is required to be small compared to the timescale 
on which the outflow rates vary. The distance traveled after N 
discrete time steps At by a parcel of gas launched at time is 



xitN, to) = ^ v(ti) At. 



(10) 



i=0 



When two parcels collide (i.e. a fast-moving parcel surpasses 
a slower moving parcel that was ejected earlier), momentum is 
conserved. It is possible that some kinetic energy is dissipated 
during a collision. Let us consider two parcels with masses mi 
and m2, pre-shock velocities vi and V2 and post-shock veloci- 
ties v[ and In order to calculate the energy loss, we perform 
a Galilean transformation to the center-of-momentum (COM) 
frame: 



V2 



vi -V = vi 

V2-V = V2, 



(11) 

(12) 



where tildes indicate quantities in the COM frame, and 

mi -\- 1712 

is the velocity of the COM. The momentum equation reads 

p = p' = 0. (14) 

We parametrize the energy loss in one collision by the factor fj^, 
defined such that 



EL = r£i 



kin? 



(15) 



where < if' < 1. Solving these equations for the velocities 
before and after the collision (assuming that no mass transfer 
occurs), and transforming back to the frame of reference of the 
observer, we have 



v[ = V±fj 



V2 = V + 7] 



m\ + m2 
^ mi(v2-vi) 



mi -h m2 



(16) 
(17) 



We choose the solution with the plus sign in Eq. (16) and minus 
sign in Eq. (17), and require the energy loss parameter fj^ to take 
values between -1 and 0. A negative 77 can thus be interpreted 
as a 'stickiness factor'. In our simulations, we have taken values 
between the two limiting cases fj = -I (the two parcels pass each 
other without any interaction) and ^ = (a fully inelastic colli- 
sion occurs). The energy loss in one collision event is expressed 
in terms of the velocity difference and fj as 



A^kin = Afikin = 



■ V2)\ 



(18) 



where ji' = mim2/(mi -h m2) is the reduced mass. 

Since a fraction if- of the center-of-mass kinetic energy is 
dissipated in every collision, the total energy loss over a longer 
period scales exponentially with the number of collisions. In 
turn, this scales linearly with the number of particles generated 



^ 77 can be interpreted as the 'coefficient of restitution' that is used in 
classical mechanics. 
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Figure 12. Results for the simulations run for HH 1042 for the stochastic {left) and periodic {right) models. The displayed stochastic result is 
selected from the 1000 generated stochastic models as it most adequately reproduces the observations. The periodic model parameters are listed 
in Tab. 6. Colors indicate the modeled density field. White contours indicate the extinction-corrected observed [Fe ii] 1643 nm flux. The result is 
displayed for increasing values of the (global) energy loss parameter, 77 = (-1, -0.5,0) from top to bottom; 77 = -1 represents a free flow while 
77 = represents a flow with fully inelastical collisions. Note the phase diff'erence of A0 = 0.65 between the red and blue lobe. 



per unit time in the simulation, which we wish to be an arbi- 
trary parameter. We describe the energy loss in terms of a global 
parameter 77: 

r^^Sigm\lf'^''-'^'\ (19) 

with 0.5 Ti the characteristic timescale for variations in the flow 
velocity and the time step. In our simulations, Ti is the longest 
period in the simulated outflow rate, Eq. (8). In this definition, if' 
is a measure of the fraction of kinetic energy dissipated over one 
characteristic timescale. Using this global parameter, the results 
of the simulation is independent of the time step used. 

5.3. Results of simulation 

The simulations, which use the input and model described in the 
previous sections, consist of 500 time generations, running from 
^tan/ = -320 yr in the past up to the present day (^tan/ = 0). 
To compare with observations, the absolute distances and ve- 
locites are converted to observable parameters {x sin / [''],v cos /) 
by adopting d = 0.7 kpc (Liseau et al., 1992). The fit parameters 
are thus degenerate with the inclination angle /. 

Because of the limitations of our simulations - the ballistic 
approach, the inclusion of only one dimension, and the lack of a 
physical model that describes the density field, energy loss and 
emergent line emission - it may only serve as an interpretative 



model. We thus limit ourselves to qualitative comparisons be- 
tween data and simulations^. 

In Fig. 12 the results of the simulations are shown for both 
the stochastic and periodic models, for three diflTerent values of 
7], increasing from 77 = -1 up to 77 = 0. Values of the parameter 
T] close to -1 lead to a larger spread in velocity, while 77 close to 
zero leads to an increased density in the knots. Both the stochas- 
tic and the periodic outflow rates can produce a shocked density 
structure in the jet. 

The periodic model leads to a better representation of the ob- 
servations than the stochastic model. The stochastic model dis- 
played in Fig. 12 is the one that out of the 1000 generated models 
most adequately represents the observations. While few of the 
stochastic models successfully reproduce the exact locations of 
the observed knots, the typical distance between them and their 
varying intensity is comparable to the observations. 

The periodic rate produces a larger spread in velocity (e.g. 
the 'saw-tooth' patterns in position- velocity space) which is gen- 
erally not generated by the stochastic model. It reproduces the jet 
emission pattern in position- velocity space reasonably well. The 
locations of knots B, C, D and E in the blue lobe (and their sub- 
structure) are reproduced. The location of knot G (see Fig. 1) is 
also correctly predicted by the simulation. 



^ We have run several tests to quantitatively compare the models to 
each other and to the data. However, no clear quantitative criteria can 
be established given the steep gradients present in the two-dimensional 
images. 
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The validity of this interpretative model and its implications 
for constraints on the launching mechanism are further discussed 
in Sees. 6.3 and 6.4. 



6. Discussion 

In this section, we discuss the validity and context of our results. 
Sect. 6.1 summarizes our findings on the mass and inclination 
angle of the two objects. In Sect. 6.2 we compare the obtained 
physical characteristics of both sources with other YSOs across 
the mass spectrum. Finally, the method described in Sect. 5 al- 
lows us to draw conclusions on the timescales, velocities, sym- 
metry and collisions in the HH 1042 jet. These are sunmiarized 
and discussed in Sect. 6.3. Some scenarios describing the nature 
of the jet-launching mechanism within the constraints obtained 
from our analysis are put forward in Sect. 6.4. 

6.1. Mass and inclination angle 

Two important parameters are not well constrained: the mass 
of the central object, and the inclination of the disk-jet system. 
Knowing the mass of the central object is relevant because op- 
tical jets are rarely found around the more massive YSOs. The 
inclination is important because it remains as a free parameter in 
the mass loss estimates, as well as in the velocities, lengths and 
timescales derived in the kinematics analysis. 

Lacking proper motion measurements, the best estimate of 
the inclination angle of the jets is represented by the inclina- 
tion of the disk (assuming that the disk rotation axis is paral- 
lel to the jet). From fitting the CO-bandhead feature at 2.3 //m, 
the inclination angle of 08576nr292 is estimated at / = 27°_\^ 
by Bik & Thi (2004), and at / = 17.8°+^/ by Wheelwright et al. 
(2010). However, these results are dependent on the adopted stel- 
lar mass, because it sets the width of the modeled line profiles. 
The inclination angle in these models, keeping all other param- 
eters fixed, then scales as sin / oc M^^. One should thus attempt 
to refine the estimate of M* to get a more accurate estimate of /. 

The value M* = 6.1 M© adopted in both studies mentioned 
above, based on the classification as an early B-type star by 
Bik et al. (2006), is probably an overestimate. It does not take 
into account the intrinsic infrared excess that is evident from 
the SED-fitting by EUerbroek et al. (201 1). From the extinction- 
corrected SED and the distance, and assuming the photosphere 
to be veiled by flat continuum emission from a disk, the stellar 
luminosity is expected to be between 10 and 70 L©, equivalent 
to an A5V-B9V star with zero-age main sequence (ZAMS) ra- 
dius. The lower limit is set by the luminosity integrated over 
the spectral range observed with X-shooter and Spitzer/IRAC 
(EUerbroek et al., in prep.); the upper limit is the brightest pho- 
tosphere which could be 'hidden' under the dereddened disk 
SED (EUerbroek et al., 2011). The eff'ective temperature prob- 
ably is lower than the ZAMS value as PMS stars have bloated 
radii (e.g. Palla & Stabler, 1993; Ochsendorf et al., 2011). At 
its accretion rate (Mace ~ 10"^ M© yr~^), the star is expected 
to arrive on the main sequence while still accreting (Yorke & 
Sonnhalter, 2002). However, the star is not expected to accrete 
more than a fraction of a solar mass over the remaining disk 
lifetime (Hartmann et al., 1998). Adopting the mass-luminosity 
relation for PMS stars (Palla & Stabler, 1993) and the consider- 
ations mentioned above, we tentatively estimate the final stellar 
mass of 08576nr292 at ~ 2 - 5 M©. Accounting for this large 
uncertainty in mass, we refine the estimate of the inclination an- 



gle to / ~ 20° - 30°, which incidentally is consistent with the 
estimate of Bik & Thi (2004). 

The spectrum of 08576nr480 does not contain signatures 
of its embedded stellar source other than the Ca ii , O i and 
Bry lines, and the CO bandhead feature. Extinction prevents a 
classification of the optical SED, while at longer wavelengths 
(3 - 10//m) the emission is probably contaminated by the sur- 
rounding cloud, preventing a classification of its infrared spec- 
trum (EUerbroek et al., in prep.). The jet appears to be quite in- 
clined with respect to the line of sight due to the visibility of 
its two lobes and the bow-shock feature. For this reason, we 
have adopted / = 60°^^^ in Sect. 4.3. The CO bandhead of 
08576nr480 is quite steep; a high inclination thus implies a low 
mass (M < 1) of the central object. The appearance of the jet 
is comparable to HH 1042 in terms of excitation conditions, and 
that more H2 and H i lines are present in the nebular spectrum, 
suggesting the source is deeply embedded in the cloud and in 
an early evolutionary stage. This is consistent with an underesti- 
mated accretion rate (see Sect. 6.2). 

From the ^-band spectrum of HH 1043, it is apparent that 
a second continuum source is present at ~ 0.8'' (600 AU) west- 
ward of 08576nr480. Since the spectrum does not contain any 
stellar features, it cannot be said whether this is a companion or 
background star or scattered light of 08576nr480 oflf a clump in 
the jet. 

6.2. Physical properties, Mjet and M^cc- 

The degree of excitation in both jet spectra is high compared 
to similar objects. Fig. 13 depicts a diagram of two line ratios: 
Ka/iS 11] and [N ii]/[0 i], both known to be tracers of excitation 
(see e.g. Hartigan et al., 1994). These line ratios in HH 1042 
and HH 1043 are larger than those observed in 'high-excitation' 
spectra (Raga et al., 1996), yet they are still within the range pre- 
dicted by plane-parallel shock models by Hartigan et al. (1987), 
given the large number of free parameters that these models de- 
pend on. However, the [S iii] lines in HH 1042 and HH 1043 are 
anomalously strong: two orders of magnitude stronger than in 
both models and observations of shocked jets. Apart from inter- 
nal shocks, an external source of radiation may also contribute 
to these high ionization conditions. Candidates for this are either 
the relatively hot central source or the nearby O stars in RCW 36, 
or both. 

Fig. 14 displays the observed mass accretion and mass out- 
flow rates of a sample of YSOs over a broad range in mass, evo- 
lutionary stage, and other properties. A correlation is seen be- 
tween Mace and Mjet, albeit with a large scatter. Most sources are 
accretion-dominated (Mjet/Macc < 1), while the accretion and 
outflow rates are highest in the more massive (HAeBe), young 
(Class O/I) and extreme (FU Ori, continuum stars) objects. The 
mass outflow rates Mjet are in some cases determined from spa- 
tially resolved jets, but mostly determined by unresolved on- 
source forbidden emission line diagnostics (cf. Sect. 4.3). The 
typical errors in the measurements of Mjet and Mace are quite 
large, due to the uncertainty in the heating mechanisms and 
structure of the shocks, the structure of the accretion region, and 
the stellar parameters. These uncertainties can account for most 
of the scatter in Fig. 14. Finally, the estimates of Mjet and Mace 
difl'er greatly depending on the literature consulted, which may 
be due to the variable nature of the accretion/ejection process 
and/or the large uncertainties in the diagnostics used. 

The accretion and outflow rates of HH 1042 (08576nr292) 
are consistent with it being more massive than a CTTS and in 
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Figure 13. Line ratios of HH sources with 'low-' (open circles), 
'intermediate-' (filled circles) and 'high-excitation' (crosses) spectra 
(data and definitions adopted from Raga et al., 1996). Red symbols rep- 
resent the brightest knots in HH 1042 (knot E) and HH 1043 (knot B). 
The shaded region denotates the ratios predicted by the plane-parallel 
shock and bow-shock models in Hartigan et al. (1987). The high exci- 
tation in HH 1042 and HH 1043 can be explained by photoionization 
and/or -excitation by an external source of radiation. 



a relatively early stage of evolution. The rates are higher than 
those observed in CTTS and are comparable with HAeBe-stars 
and the continuum stars, to which 08576nr292 bears many spec- 
troscopic similarities. The accretion luminosity is of the same 
order as the stellar luminosity, which is consistent with mod- 
els for intermediate-mass stars with high accretion rates (Palla 
& Stabler, 1993). Our results imply a Mjet/Macc ratio of ~ 0.1, 
consistent with values predicted by magneto-centrifugal forces 
(see e.g. Cabrit, 2009). 

For HH 1043 (08576nr480), the estimated ratio M-^JM^cc 
exceeds unity. This indicates either an 'outflow-dominated' sys- 
tem, or an underestimate of the accretion rate due to the accretion 
region being too embedded or obscured to be probed by the di- 
agnostics used in Sect. 4.4 (see e.g. the discussion in Bacciotti 
et al., 201 1). Based on the previous, 08576nr480 is likely a Class 
O/I object. 

6.3. Validity of the kinematics analysis and Interpretative 
model 

The results of Sect. 5 allow us to derive quantitative constraints 
on the launching mechanism of the HH 1042 jet (and to a lesser 
extent the HH 1043 jet), and also to comment qualitatively on its 
general outflow properties. First of all, the dynamical timescale 
(to be interpreted as the minimum jet lifetime, based on the mean 
velocity and spatial extent of the jets, taking into account the 
uncertainty in inclination) is 500 - 1000 yr for the blue lobe of 
HH 1042. The red lobe is visible up to 300 - 600 yr. For HH 1043 
the dynamical timescales are 100 - 300 yr (blue lobe) and 100 - 
400 yr (red lobe). The periods associated with the spatial in- 
tervals between the knots and the modulations of the outflow 
velocity are in the order of 10 - 100 yr in both jets. More specif- 



Figure 14. Observed mass outflow rate Mjet versus accretion rate M^cc 
of HH 1042 and HH 1043 (red squares) compared with difi'erent classes 
of YSOs associated with jets and outflows. The M^cc value of some 
HAeBe stars (grey triangles) was determined from the Ra luminosity 
in Hillenbrand et al. (1992). The orange symbols indicate 'continuum 
stars', objects classified as CTTS (orange circles) or HAeBe stars (or- 
ange triangles) with a high degree of veiling whose spectral type is very 
uncertain (see Calvet & Gullbring, 1998; Hernandez et al., 2004); these 
objects are spectroscopically similar to 08576nr292, the driving source 
of HH 1042. The dashed lines indicate the Mjet/M^cc = 0.01 and 0.1 
range, which is predicted by jet formation models. 
References: Class O/I objects: Hartigan et al. (1994); Bacciotti & 
Eisloffel (1999); Podio et al. (2006, 2011, 2012); Antoniucci et al. 
(2008); Cabrit (2009); FU Orionis objects: Hartmann & Calvet (1995); 
Calvet (1998); HAeBe stars: Levreault (1984); Boehm & Catala (1994); 
Nisini et al. (1995); Garcia Lopez et al. (2006); Wassell et al. (2006); 
Melnikov et al. (2008); Mendigutia et al. (2011); Donehew & Brittain 
(2011); Liu et al. (2011); CTTS: Hartigan et al. (1995); Gullbring et al. 
(1998); Muzerolle et al. (1998a); Coffey et al. (2008); Cabrit (2009); 
Melnikov et al. (2009); Brown dwarfs: Whelan et al. (2007); Herczeg 
& Hillenbrand (2008); Bacciotti et al. (2011). 

ically, the shape of HH 1042 in position- velocity space can be 
reconstructed reasonably well with periodic modulations of the 
outflow velocity v(t) (with period Ti tan / = 83 yr) and the mass 
outflow rate m(t) (7^2 tan/ =12 yr). The stochastic model pro- 
duces similar knotted structure but a smaller velocity spread. The 
faint spike in the flux periodogram at Ttan/ ~ 80 yr (Fig. 11) 
may reflect a periodic signal in the outflow rate. Incidentally, it 
coincides with the best-fit value Ti tan / of the longest period in 
the periodic model. 

Instead of modulating m(t), one could also reproduce the 
knotted structure by superposing a second small amplitude mod- 
ulation on v(t), as done in e.g. Raga et al. (2012). The resolu- 
tion of our data does not allow us to distinguish between either 
a single-mode velocity and mass outflow modulation or a two- 
mode velocity modulation. In Fig. 9 (top) it can be seen that at 
the location of the knots, the velocity drops and the width of the 
individual components in the spectral profile increases. This sug- 
gests that they consist of diflferent components scattered in veloc- 
ity, which is not an uncommon observation in jets (Hartigan & 
Hillenbrand, 2009). While the detection of multiple shock fronts 
would put constraints on the launching mechanism, the spectral 



17 



L.E. EUerbroek et al.: The history of HH 1042 and HH 1043 



and spatial resolution of our observations are insufficient to test 
for their existence. 

The mean velocity in HH 1042 is comparable in the blue 
and red lobe, f;o cos / ~ 170 km s"^ The amplitude of the mod- 
ulations on this velocity is ficos/ ~ 40 km s"^ . HH 1043 
has mean velocities that are somewhat more asymmetric (80 
and 60 km s"^ in the blue and red lobes, respectively), but it 
could be that the system is moving with 10 km s"^ toward us 
with respect to the ambient cloud. Both the velocities and the 
timescales are similar to those derived for other jet sources (e.g., 
HH 30, HH 34, HH 110, HH 111, HH 444, Hen 3-1475; Raga 
et al., 2001, 2002a,b, 2004, 2010, 2012; Velazquez et al., 2004; 
Esquivel et al., 2007). Moreover, in most of these sources, like in 
the case of HH 1042, multiple modes with different timescales 
of order 10-1000 yr are found. Two timescales found in a sin- 
gle source typically have a ratio of order 10. These results may 
indicate that the mechanism responsible for jet launching is in- 
deed intrinsically variable, and similar mechanisms causing vari- 
ability may be at work in all HH objects across the stellar mass 
range. 

In both HH 1042 and HH 1043, the knots in the blue lobe 
do not have obvious counterparts in the red lobe. Consequently, 
there appears to be no correlation between ejection events in the 
blue and red lobes. A phase offset of A0 = 0.65 was adopted in 
the periodic models, but the parameter is not well constrained; 
a value of 0.5 cannot be excluded. The latter value would be a 
natural consequence of a precession movement of the jet (see the 
next subsection). 

The simple ballistic model presented in this paper can quali- 
tatively reproduce the structure observed in the jets. Shocks form 
along the jet wherever fast material catches up and collides with 
slower material. The local increase in density then results in 
an increased line emissivity (e.g. Bacciotti & Eisloffel, 1999). 
However, particularly in and beyond shock fronts (knot E), there 
is a mismatch in both velocity and flux level, indicating that a 
more sophisticated treatment of collision physics is needed. It 
is clear that some fraction of the kinetic energy is dissipated in 
collisions and either ionizes or excites the material, which heats 
and/or re-emits in different emission lines (Hartigan et al. e.g. 
1994; see Fig. 5). The inclusion of a gas-dynamical model pre- 
dicting line emissivities and fitting the data to more than one line 
species would likely improve these simulations, but is beyond 
the scope of this paper. 

Models (both stochastic and periodic) with the energy loss 
parameter 77 close to -1 barely have collisions and are therefore 
better suited to explain the velocity spread. However, it is evi- 
dent from the variation in line emission (Eig. 5) that ionization 
conditions vary along the jet, probably caused by shocks, which 
suggests that collisions do occur in the flow. In reality, the pa- 
rameter 7] may not have a single, global value as assumed in our 
models, but rather be dependent on local physical variables like 
density, temperature and degree of ionization. 

6.4. The nature of the jet-launching mechanism 

The outcome of the simulations in Sect. 5 suggest that while 
the HH 1042 outflow rate may not be purely periodical, a quasi- 
periodic variation on the quoted typical timescales cannot be ex- 
cluded. A straightforward explanation of the variability of veloc- 
ities in the jet is an intrinsic variability of the accretion/ejection 
mechanism at the base of the jet. Variability of a quasi-periodic 
nature (i.e. with a characteristic, but not fixed timescale) can also 
be related to the piling up and subsequent release of material in 



the inner disk, or to the reconfiguration of the stellar magnetic 
field due to the rotation of the disk with respect to the star. 

In either of these scenarios, the origin of the phase shift is 
still unclear. It may be explained by a launching mechanism that 
does not produce knots in the red and blue lobes simultaneously. 
An alternative explanation for variability of the jet velocity is 
precession of the jet, which induces a half-period phase shift be- 
tween the blue and red lobes. Erom the collimation seen in the 
SINEONI data (Eig. 1) the precession angle is constrained to be 
less than 10°. The observed velocity mode vi could then be ex- 
plained if the inclination of the system is ~ 60°, corresponding to 
a < 1 M© star. Considering that the object is likely more massive, 
this scenario is not favored. A similar explanation for the phase 
lag is the movement of the jet on a binary orbit, which is not per- 
fectly perpendicular to the jet axis. However, the values for the 
period and amplitude of this orbit render an unrealisticly high 
companion mass. Additional explanations for asymmetries in the 
launching mechanism are different ISM conditions or magnetic 
field configurations (Matsakos et al., 2012). However, these gen- 
erally lead to a systematic asynmietry in the velocity of the blue 
and red lobes, while we observe a temporal asymmetry in the 
velocity modulations. 

7. Summary 

We have presented a detailed study of the physical properties 
and kinematics of two newly discovered jets in RCW 36. These 
are interesting objects because of the visibility of both of their 
lobes, their being located in a massive star-forming region, the 
asymmetry in their kinematic structure, and the central source 
of HH 1042 being an intermediate-mass YSO. Below, the most 
important findings of this study are sunmiarized. 

- Both HH 1042 and HH 1043 display a shocked kinematic 
structure, with high ionization conditions in their terminal 
knots and at shock fronts along the jet (knot E in HH 1042). 
HH 1042 shows a 'saw-tooth' pattern in position- velocity 
space, while the knots of HH 1043 are more widely sepa- 
rated spatially. The HH 1042 red lobe disappears into or be- 
yond the molecular cloud, while its blue lobe terminates in 
an interaction region with the ISM (knot G). The HH 1043 
blue lobe terminates in a clear bow- shock shaped structure 
(knot B). 

- The electron density in both jets is ~ 10^-10^ cm~^. The 
electron temperature is ~ 10"^ K close to the source, then 
decreases further away from the source, similarly to what is 
found for low-mass YSOs. The ionization fraction is esti- 
mated to be low on-source (< 0.025 in HH 1042 A), while it 
is very high in the outer knots (~ 0.7 in HH 1042 E). 

- The on-source extinction of 08576nr292 measured in the jet 
is several orders of magnitude lower than the value mea- 
sured in the accretion column, suggesting the presence of cir- 
cumstellar dust close to the star. The extinction in HH 1042 
steeply increases in the red lobe, consistent with its disap- 
pearance into or behind a molecular cloud. 

- Both jets have a mass outflow rate of order Mjet ~ 
10"^ Mo yr~\ comparable to what is found for HH objects 
from low-mass stars and more evolved intermediate-mass 
stars (HAeBe). The estimated accretion rate of HH 1042 
(Mace = l-10!o 2i ^ yr"^) is high compared to 
CTTS accretion rates, consistent with a higher mass and/or 
the fact that its photosphere is heavily veiled by contin- 
uum emission from an optically thick accretion disk. The 
accretion rate of 08576nr480 (HH 1043) is more modest 
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(Mace = 5.50+1^^ X 10"^ Mo yr-i), yet it is likely under- 
estimated as the source is embedded. 

- For HH 1042, the combined optical and near-infrared diag- 
nostics yield Mjet/Macc ~ 0.1, which is consistent with val- 
ues predicted by magneto-centrifugal forces (Cabrit, 2009), 
and with values observed in other, comparable sources. 

- The kinematic structure of HH 1042 can be well simulated 
by an interpretative model with as inputs a stochastic or pe- 
riodically variable outflow rate and a ballistic flow. The ve- 
locities (100 - 200 km s"^) and derived timescales (10-100 
yr) are typical of those seen in other HH objects. 

- A (quasi-)periodic signal is detected in the data at T tan i ~ 
80 yr (i.e. T ~ 100 - 200 yr). The stochastic and periodic 
models both reproduce the knotted structure and apparent pe- 
riodicity. The periodic model better reproduces the velocity 
spread in the data. 

- The velocity spread is best reproduced by models with 
T] = -I (i.e. a small amount of kinetic energy dissipation). 
However, the amount of energy loss likely depends on local 
conditions in the flow. 

- The asymmetry between the positions and apparent launch- 
ing times of the knots in the red and blue lobe of HH 1042 is 
well fitted by a phase shift of the outflow rate of 0.65 period, 
while a phase shift of 0.5 period cannot be excluded. Possible 
explanations for such a phase shift are an intrinsically asym- 
metric launching mechanism, or a precession movement of 
the jet. 

- 08576nr292, the driving source of HH 1042, is an 
intermediate-mass YSO based on its jet and accretion prop- 
erties. It resembles the class of 'continuum stars': HAeBe- 
like stars with a photosphere veiled by extreme accretion ac- 
tivity. 

We have put forward various scenarios that might lead to 
the observed structure in the HH 1042 jet. Further high angular 
resolution and/or proper motion data may help to distinghuish 
between these. 
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L.E. EUerbroek et al.: The history of HH 1042 and HH 1043 
Table A.l. Integrated fluxes per knot in HH 1042 (blue lobe). 



Flux (10-^^ erg s"^ cm"^) 



Line 


A (nm) 


A 




B 






C 






D 






E 






F 






[0 I] 


630.0 


9.28 


+ 0.11 


1.43 


+ 


0.09 














2.13 


+ 


0.10 








[0 I] 


636.4 


2.21 


+ 0.10 
































[Nn] 


654.8 












2.61 


+ 


0.06 


1.41 


+ 


0.04 


4.83 


+ 


0.07 


3.37 


+ 


0.07 


Ha 


656.3 


16.9 


+ 0.1 


5.54 


+ 


0.07 


19.8 


+ 


0.1 


10.9 


+ 


0.1 


36.7 


+ 


0.1 


23.3 


+ 


0.1 


[Nil] 


658.3 


0.50 


+ 0.05 


2.07 


+ 


0.04 


7.99 


+ 


0.05 


4.45 


+ 


0.04 


14.4 


+ 


0.1 


7.26 


+ 


0.07 


[S II] 


671.6 


0.52 


+ 0.05 


0.95 


+ 


0.04 


1.40 


+ 


0.04 


0.79 


+ 


0.03 


2.92 


+ 


0.05 


1.48 


+ 


0.05 


[S II] 


673.1 


1.15 


+ 0.04 


1.44 


+ 


0.04 


2.48 


+ 


0.04 


1.43 


+ 


0.03 


5.43 


+ 


0.05 


2.88 


+ 


0.06 


[Fe II] 


715.5 


4.35 


+ 0.04 


0.74 


+ 


0.03 


0.76 


+ 


0.03 


0.47 


+ 


0.02 


1.56 


+ 


0.03 


0.73 


+ 


0.03 


[Ca II] 


729.1 


1.77 


+ 0.05 


0.31 


+ 


0.04 


0.56 


+ 


0.04 


0.35 


+ 


0.03 


1.09 


+ 


0.04 


0.46 


+ 


0.04 


[0 II] 


732.0 












0.79 


+ 


0.04 


0.36 


+ 


0.03 


1.33 


+ 


0.04 


0.84 


+ 


0.04 


[Ni II] 


737.8 


1.66 


+ 0.04 


0.52 


+ 


0.03 


0.80 


+ 


0.03 


0.50 


+ 


0.02 


1.73 


+ 


0.03 


0.74 


+ 


0.03 


[Fen] 


745.3 


2.07 


+ 0.04 


0.31 


+ 


0.03 














0.58 


+ 


0.02 








Oi 


844.6 


5.69 


+ 0.05 


1.44 


+ 


0.05 














1.60 


+ 


0.02 








Pa- 14 


859.8 


1.00 


+ 0.04 
































[Fell] 


861.7 


10.3 




2.20 


+ 


0.04 


2.28 


+ 


0.03 


1.48 


+ 


0.02 


5.75 


+ 


0.03 


2.49 


+ 


0.03 


Pa-12 


875.0 
























0.73 


+ 


0.03 








Pa- 11 


886.3 


2.98 


+ 0.06 


0.34 


+ 


0.05 














1.01 


+ 


0.03 


0.81 


+ 


0.03 


[Fe n] 


889.2 


1.94 


+ 0.05 


0.69 


+ 


0.05 


0.90 


+ 


0.03 


0.52 


+ 


0.02 


2.16 


+ 


0.03 


1.21 


+ 


0.03 


Pa-10 


901.5 


5.46 


+ 0.07 


0.98 


+ 


0.06 


0.97 


+ 


0.03 


0.54 


+ 


0.02 


2.00 


+ 


0.03 


1.65 


+ 


0.03 


[Fe n] 


905.2 
























2.62 


+ 


0.03 








[S III] 


906.9 


1.28 


+ 0.07 


1.24 


+ 


0.07 


6.75 


+ 


0.04 


4.47 


+ 


0.03 


15.3 






17.6 


+ 


0.1 


Pa-9 


922.9 


5.02 


+ 0.10 


1.10 


+ 


0.09 


1.74 


+ 


0.05 


0.93 


+ 


0.03 


4.06 


+ 


0.05 


3.22 


+ 


0.05 


[S III] 


953.1 






6.90 


+ 


0.15 


40.8 


+ 


0.1 


24.8 


+ 


0.1 


88.9 


+ 


0.2 


102. 






Pa-8 


954.6 


16.0 


+ 0.2 




















9.34 


+ 


0.09 








Pa^ 


1004.9 


12.1 


+ 0.2 




















9.36 


+ 


0.12 








Hex 


1083.0 


3.50 


+ 1.29 








7.16 


+ 


1.79 


9.72 


+ 


1.34 


54.7 


+ 


2.5 


92.0 


+ 


14.6 


Pay 


1093.8 


32.1 


+ 1.5 


3.17 


+ 


1.20 


9.41 


+ 


1.51 


5.80 


+ 


1.12 


20.8 


+ 


2.2 


23.7 


+ 


18.4 


[Feii] 


1256.7 


57.5 


+ 1.0 


31.7 


+ 


0.8 


59.5 


+ 


1.0 


45.5 


+ 


0.8 


156. 


+ 


1. 


83.6 


+ 


11.8 


PayS 


1281.8 


81.0 


+ 1.5 


10.1 


+ 


1.0 


27.5 


+ 


1.3 


16.7 


+ 


1.0 


67.1 


+ 


1.5 


83.6 


+ 


12.1 


[Fe n] 


1294.3 


14.7 


+ 1.6 


4.03 


+ 


0.83 


9.00 


+ 


0.97 


6.39 


+ 


0.74 


22.0 


+ 


1.1 


13.6 


+ 


11.1 


[Feii] 


1320.6 


47.9 


+ 1.4 


13.7 


+ 


1.1 


24.9 


+ 


1.3 


18.9 


+ 


1.0 


64.7 


+ 


1.5 


37.0 


+ 


22.8 


[Fen] 


1533.5 


12.1 


+ 0.8 


2.38 


+ 


0.47 


7.58 


+ 


0.54 


6.61 


+ 


0.40 


21.4 


+ 


0.6 








[Fe n] 


1599.5 


18.8 


+ 0.8 


4.04 


+ 


0.53 


7.74 


+ 


0.70 


5.81 


+ 


0.79 


22.1 


+ 


0.8 


15.9 


+ 


9.0 


Br- 13 


1610.9 


15.7 


+ 0.7 
































Br-12 


1640.7 


18.0 


+ 0.6 
































[Fell] 


1643.5 


59.7 


+ 0.5 


39.3 


+ 


0.4 


75.0 


+ 


0.5 


57.9 


+ 


0.4 


205. 


+ 


1. 


126. 


+ 


7. 


[Fen] 


1663.8 


5.90 


+ 0.64 


2.19 


+ 


0.39 


4.56 


+ 


0.43 


3.42 


+ 


0.33 


12.8 


+ 


0.5 


9.14 


+ 


8.15 


Br- 11 


1680.7 


26.1 


+ 0.7 




















4.44 


+ 


0.52 








Br- 10 


1736.2 


52.3 


+ 0.7 
































Br-9 


1817.4 


183. 


+ 14. 


17.1 


+ 


7.9 


21.3 


+ 


9.3 


7.92 


+ 


6.87 


30.3 


+ 


10.2 








Br-8 


1944.6 


594. 


+ 3. 


23.0 


+ 


2.3 


20.2 


+ 


2.8 


10.4 


+ 


2.1 


42.7 


+ 


3.1 


72.1 


+ 


51.9 


Hex 


2058.1 


15.6 


+ 1.5 




















6.83 


+ 


0.92 








Bry 


2165.5 


103. 


+ 1. 


6.45 


+ 


0.36 


13.2 


+ 


0.4 


7.57 


+ 


0.27 


31.2 


+ 


0.4 


49.5 


+ 


9.6 



L.E. EUerbroek et al.: The history of HH 1042 and HH 1043 
Table A.2. Integrated fluxes per knot in HH 1042 (red lobe). 



Line 


A (nmj 


A 


t71tt^ ^^ A- 
rlUX (iU 

rJ 


erg s ^ cm ^) 

r^t T\t T2' 

L. D L 


[0 i] 


630.0 


r\ r\i—i , r\ i r\ 

2.07 + 0.10 






[0 i] 


636.4 








[N n] 


654.8 








Ha 


656.3 


107. 


2.26 + 0.05 




[N ii] 


658.3 








[S II] 


671.6 


f\ no 1 f\ A/1 

0.78 + 0.04 


1 AC 1 A AC 

1.05 + 0.05 




[S II] 


673.1 


0.88 + 0.04 


1 /I O 1 A A /I 

1.43 + 0.04 




[Fe ii] 


715.5 


0.77 + 0.03 






[Ca ii] 


729.1 


1 An 1 f\ f\ A 

\A1 ± 0.04 






[0 II] 


732.0 


^ OO 1 f\ AO 

2.88 + 0.03 






[Ni ii] 


737.8 








[Fe ii] 


'7/1 C O 

745.3 








O I 


844.6 


12.1 


OO 1 A AO 

2.33 + 0.03 




ra-14 


ocn o 


O 1 O 1 A AO 

3.18 + U.Uz 






[Fe II] 


861.7 


OO 1 A A'^ 

2.83 + 0.02 


A AO 1 A A'^ 

0.98 + 0.02 




Pa- 12 


875.0 


/C AT 1 A AO 

6.07 + 0.03 






Pa- 11 


886.3 


1 1 O 

11.8 






[Fe ii] 


889.2 








Pa- 10 


901.5 


11.0 






[Fe ii] 


905.2 








[S III] 


906.9 








Pa-9 


922.9 


14.6 






[S m] 


953.1 








Pa- 8 


954.6 


51.2 + 0.1 






Pa^ 


1004.9 


39.5 + 0.1 






He I 


1083.0 


62.1 + 1.5 






Pay 


1093.8 


59.8 + 1.1 






[Fe II] 


1256.7 


19.2 + 1.0 


30.5 + 1.0 


o OO 1 A on 

8.82 + 0.87 


Pa/3 


1281.8 


166. + 1. 


O CO 1 A AA 

3.68 + 0.90 




[Fe ii] 


1294.3 


A 'T7 1 A O/C 

4.27 + 0.86 


O OO 1 A /CO 

3.33 + 0.62 




[Fe ii] 


1320.6 


19.4 + 1.0 


10.7 + 0.7 


O A/1 1 A cn 

3.04 + 0.57 


[Fe II] 


1533.5 


24.1 + 0.7 


^ 1 A 1 A OO 

6.19 + 0.32 


1 A 1 1 A O C 

1.91 + 0.25 


[Fe ii] 


1599.5 


6.16 + 0.59 


O OT 1 A O C 

3.37 + 0.35 




Br- 13 


1610.9 


O O ^ 1 A /I A 

3.36 + 0.49 






Br- 12 


1640.7 


A TA 1 A /I A 

9.79 + 0.49 






[Fe II] 


1643.5 


21.9 + 0.5 


40.7 + 0.3 


lOOiAO OTI lAOO 

13.3 + 0.2 8.71+0.28 


fT^f^ ttI 

L^e iij 


IOOj.O 








Br- 11 


1680.7 


9.83 + 0.53 






Br- 10 


1736.2 


21.7 + 0.5 






Br-9 


1817.4 


2797. + 268. 


268. + 129. 




Br-8 


1944.6 


268. + 2. 


5.52 + 0.92 




He I 


2058.1 


13.9 + 0.9 






Bry 


2165.5 


99.8 + 0.6 
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L.E. EUerbroek et al.: The history of HH 1042 and HH 1043 
Table A.3. Integrated fluxes per knot in HH 1043 (blue and red lobes). 



Flux (10-^^ erg s"^ cm-^) 
Blue lobe Red lobe 



Line 


X (nm) 


A 


B 


A' 


B' 


C 








[0 I] 


630.0 


















[0 I] 


636.4 


















[Nn] 


654.8 


















Hof 


656.3 


30.4 + 0.1 


40.3 + 0.1 


25.2 + 0.1 


21.9 + 0.1 


6.54 + 0.06 


3.83 


+ 


0.06 


[Nil] 


658.3 


6.97 + 0.07 


14.3 + 0.1 


9.20 + 0.07 


7.09 + 0.07 


1.60 + 0.04 


1.12 


+ 


0.04 


[S n] 


671.6 




2.31 + 0.05 


1.10 + 0.05 


0.79 + 0.05 










[S 11] 


673.1 


1.50 + 0.05 


4.47 + 0.05 


2.05 + 0.05 


1.38 + 0.05 


0.29 + 0.03 








[Fe II] 


715.5 


















[Ca II] 


729.1 


















[0 II] 


732.0 


1.75 + 0.05 


3.56 + 0.05 


2.31+0.03 


1.27 + 0.03 


0.28 + 0.02 








[Ni II] 


737.8 


















[Feii] 


745.3 


















Oi 


844.6 


















Pa- 14 


859.8 


















[Feii] 


861.7 


1.88 + 0.03 


2.09 + 0.02 


1.93 + 0.03 


1.18 + 0.03 










Pa-12 


875.0 


















Pa- 11 


886.3 


















[Feii] 


889.2 


















Pa- 10 


901.5 




2.07 + 0.03 


1.76 + 0.03 


1.56 + 0.03 










[Fe II] 


905.2 


















[S III] 


906.9 


19.3 + 0.1 


19.6 + 0.1 


33.6 + 0.1 


40.7 + 0.1 


14.2 + 0.0 


8.20 


+ 


0.04 


Pa-9 


922.9 


3.77 + 0.07 


4.15 + 0.06 


3.04 + 0.06 


2.64 + 0.05 










[S III] 


953.1 


138. 


169. 


109. 


101. 


37.1+0.1 


24.9 


+ 


0.1 


Pa-8 


954.6 








8.26 + 0.09 










Pa^ 


1004.9 






11.7 + 0.1 


11.0 + 0.1 


3.44 + 0.07 








Hex 


1083.0 


198. + 2. 


212. + 2. 


172. + 2. 


155. + 2. 


66.9 + 1.0 


54.5 


+ 


1.1 


Pay 


1093.8 


28.0+1.6 


27.3+1.4 


22.0+1.0 


21.2 + 0.9 


5.18 + 0.57 


4.90 


+ 


0.61 


[Fen] 


1256.7 


59.1 + 1.0 


68.9 + 0.9 


64.7 + 1.4 


55.5 + 1.2 


6.53 + 0.51 


18.2 


+ 


0.6 


Pa/? 


1281.8 


94.5 + 1.5 


84.1 + 1.4 


92.6 + 1.3 


88.0+1.0 


21.8 + 0.6 


20.2 


+ 


0.9 


[Feii] 


1294.3 


9.28 + 0.98 


10.2 + 0.9 


12.8 + 0.7 


9.11 + 0.61 




3.30 


+ 


0.55 


[Fen] 


1320.6 


24.0+1.0 


24.0+1.0 


20.9 + 0.9 


18.5 + 1.0 


2.12 + 0.45 


6.25 


+ 


0.65 


[Fe n] 


1533.5 


15.1+0.6 


16.8 + 0.6 


24.2 + 0.4 


19.6 + 0.3 


2.86 + 0.19 


7.07 


+ 


0.39 


[Fen] 


1599.5 


14.3 + 0.9 


12.4 + 0.7 


16.1 + 0.5 


11.9 + 0.3 


1.26 + 0.21 


4.75 


+ 


0.34 


Br-13 


1610.9 


5.07 + 0.47 




3.85 + 0.53 


4.09 + 0.29 




1.10 


+ 


0.27 


Br- 12 


1640.7 


5.80 + 0.38 


4.75 + 0.34 


4.52 + 0.32 


4.91 + 0.27 


1.12 + 0.17 


1.12 


+ 


0.21 


[Fen] 


1643.5 


95.4 + 0.9 


103. 


103. 


90.1 + 0.7 


10.1 + 0.2 


33.7 


+ 


0.3 


[Fen] 


1663.8 


7.13 + 0.42 


7.13 + 0.37 


9.25 + 0.41 


7.18 + 0.28 




2.95 


+ 


0.31 


Br- 11 


1680.7 


13.2 + 0.5 


12.9 + 0.4 


4.00 + 0.38 


5.20 + 0.28 


1.27 + 0.17 


1.62 


+ 


0.36 


Br-10 


1736.2 


14.4 + 0.5 


9.88 + 0.68 


12.0 + 0.3 


11.8 + 0.2 


2.91 + 0.14 


3.26 


+ 


0.27 


Br-9 


1817.4 


46.9+14.1 


35.9 + 9.3 


244. + 222. 












Br-8 


1944.6 


337. + 24. 


410. + 22. 


42.9 + 1.1 


52.3 + 0.9 


11.6 + 0.6 


13.5 


+ 


1.0 


He I 


2058.1 


47.3 + 0.9 


25.6+1.0 


31.8 + 0.9 


33.7 + 0.5 


11.5 + 0.3 


11.1 


+ 


0.6 


Bry 


2165.5 


80.9 + 0.6 


63.2 + 0.5 


74.6 + 0.4 


78.2 + 0.4 


16.5 + 0.3 


20.5 


+ 


0.4 



